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On the Scattering of Waves by a Disk 


By Avsert E. Herns, Ann Arbor, Mich., Urs:A54) 
and Rrcwarp C, MacCamy, Pittsburgh, Pa., U.S.A.®) 


1. Introduction 


This paper is concerned with the scattering of sound waves by a circular 
_ disk of so-called ‘soft’ material. The appropriate boundary value problem is 
_ the determination of a function U(r, 0, z) satisfying the equation 


1 1 
U,, + — U, +=4 Uy + U,, +B U=0 (1.1) 
_ with the condition 
4 U(r, 0,0) forr<1. (1.2) 
4 addition U(r, 6, z) is to be made up of an incident field U,(r, 6, z), which 
is prescribed, and a term U*(r, 0, z) satisfying a radiation condition at infinity. 
_. An explicit solution of this problem is possible by means of spheroidal 


harmonics [1]’) but the resulting formulae are exceedingly complicated. An 
alternative method is to formulate the problem as an integral equation. Using 


122 


Ur, 6,2)=W.9.)+4- | [Hed qe odbdo (1.3) 


with R? = 7? + 9? — 2rcos(6 — ¢) + # and where f(o, ¢) denotes the discon- 
tinuity of U, on the disk, z = 0, y <1. From (1.2) and (1.3) one obtains the 
integral equation, 

1 22 


1 7 gikR 
0= U(r, 0,0) + 45 | | Heo) RQ ededd, P= 0 ot 
a0: 


Unfortunately a solution of (1.4) is not known and one is led to the question 
of approximations. 
1) The University of Michigan, Department of Mathematics 


2) The Carnegie Institute of Technology, Department of Mathematics. 
3) Numbers in brackets refer to References, page 264. 
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When # is set equal to zero the above problem becomes that of an elec- 
trified disk in electrostatics, simple solutions of which have long been known. 
It was observed by Copson [2] that in the case k = 0 the equation (1.4) itself 
can be solved explicitly. Copson’s procedure calls for the expansion of all 
quantities in a Fourier series of the angle § and then carrying out the angular 
integrations in (1.4) to obtain simpler integral equations which can be solved. 
For the case of axially symmetric problems, that is, when U does not depend 
on 0 the method of Corson was extend by Jones [3]*) to equation (1.4) for 
arbitrary k. In this case an explicit solution is no longer possible but rather 
Jones was able to derive from (1.4) a regular Fredholm integral equation of 
second kind which is well adapted to an iteration scheme for small k. 

It is rather an involved procedure to carry out the angular integration in 
(1.4) for k + 0. It proved necessary to use some rather complicated relations — 
involving Bessel functions. Curiously enough, these complications vanished at 
the end and the resulting equation involved only elementary functions. The 
authors were led to believe that there should be also a simpler method for 
deriving this equation and this, in fact, proves to be true. The central feature 
of the problem appears to be the fact that an axially-symmetric solution of 
(1.1) is determined uniquely by its values on the axis of symmetry, 7 = 0. Such 
a determination was observed by HEnRIcI [4] and further investigated by the 
authors in [5]. In this paper we show that an integral equation analogous to 
that of Jones may be derived without the intermediate step of rewriting the 
@ integration in (1.4). Further, the analysis is simple enough to be easily 
carried over to the case of non-symmetry about 7 = 0. 

In section 2 we indicate how the general situation without symmetry can 
be reduced to a series of problems similar to the axially symmetric case. In 
section 3 the necessary material from [4] is presented. In the next two sections _ 
we outline our method of solution and in the last section we carry out some 
simple calculations for the case in which U,(r, 0, 2) represents a plane wave of 
arbitrary angle of incidence. We include an approximate formula for the 
discontinuity of the z derivative of U on the disk, as a function of angle of ' 
incidence, when the wave number 2 is small. 


2. Preliminary Remarks 


Solutions of (1.1) which are of class C) ina region are known to be analytic : 
functions of z and x = r cos 0, y = rsin 0. The particular function U satisfying : 
(1.2) will accordingly be analytic except on the ring y = 1, z =0, so that the: 
apparent singularity at y = 0 does not in reality appear. In this connection we : 
shall have frequent use of the following lemma: 


4) Jones studied the problem with the condition 0 U/Oz = 0, z = 0 but the changes to the present 1 
problem are not difficult. 
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Lemma 1. Let G(x, y, 2) be an analytic function for x? +- y? < R*. Then with 
x=~rcos0, y=rsin6o, 


co 


Gs o [A,(7, 2) cosn 0 + B,(r, 2) sinn 6], 


n=0 
where 


A,r, 2) =r" a,(r, 2), B,(r,2) =r" b,(r, 2), 
a, and b, being continuous at ry = 0. 


That G can be expanded in a Fourier series is immediate. Further since G 
is analytic, 


0 4© © 
a, 25 Sig RLY Be): OSs a Qj mn yi Cos! sin? G27, 


7=0 m=0 n=0 7=0 m=0 n=0 


Now the Fourier series for G can be obtained by expanding cos/6 and sin™6 
_ in multiple angles. However cos/@ sin™@ will contribute to cosm @ or sinn 6 
_ only if 7 + m =n. It follows that the coefficients of cos” 6 and sinn 6 involve 
. only terms in 7” or higher powers of r. 

On the basis of lemma 1 we can write for the function U, 


n=0 


er, 0,2). = ye [U(r, z) cosn0+V%(r,z)sinn 6]; U*=ru4,V°=rv", | 


(2.1) 
rv, o—. > |U, (7,2) cosn 0+ Vor, z)sinn 0]; Ul =u, Vo =r" a | 
n=0 
_ We write also, 
fe, $) = D} ¢n(g) cosm d + s,(0) sinn 6 (2.2) 


n=0 


Since the function U is to be a solution of (1.1) the functions U" and V” are 
solutions of the -. 


W,, ++ W,+W.,+ (#- 2) w= 0 23) 


2n+1 
Y 


Wp + w,+w,,+tew=0. (2.4), 
If the series (2.1) and (2.2) are substituted in the Helmholtz formula (1.4) one 
finds easily that the terms corresponding to the various Fourier coefficients 


7 
4 
4 
; 
, 
$ 
F 
: 
7 
: from which it follows that the functions w” and v” are solutions of 
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separate. In fact if we introduce the operators 


M2a—aq | /s (0) ecosn p “5, dd do, 


1 22 


N22) = =, [8 ) osinn  *— dd do, 


n 


with R? = 7 + 0? — 27 @cosd + 2, formula (1.3) becomes, 
U(r, 2) = Un(r, 2) + M,(r, 23 ¢,), V(r, 2) = Vor, z) + ML, 2; 8,). (2.5) 


We now proceed to derive one of the two key formulae in our work. We 
prove 


on 2 N+1 piku 
n\ u"(0, 2) =n! u3(0, 2) + (—! Gey | ale) —"— do, »#0, 
0 


where uw? = 9? + 2%. The same relation exists between v", vj and s,. In order 
to verify (Z) we start with the observation that U” = 7” 4” implies, : 


on UN! OUR 
nl ut(0, 2) = S08 


n! u*(Q0, 2) = ——— ; . 
( ) Or” |r=0 Oy \r=0 


Thus we are required to compute first the successive derivatives with respect 
to y of M,(r, z; c,). In order to carry out these differentiations we write, 


—- = F(B+2), p=r+0?—2rocosd. 


We have then 


B,=2r—Zecosd, B,,=2 and eT =0 for m = 3% 
Then we assert, 
as gikR [n/2] 4 
4 FUT Bes?) Gl mee] (2.6). 


We prove (2.6) by the method of induction. It is clearly true for n= 1. Let us 
assume now that it is true for any m > 1. Then 


anti e ikR [n/2] [n/2] 


soe ae “ee Cf Peat) Bn-2)-+1 8) ee 3,4 cr Fat peste fi (n — 


since 0° B/dr3 = 0. Each term in the second sum has the same form as a term — 
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in the first sum except the last one. This one is 
Gia Fo ee Bee — 2 [n/2)) , 


Now if 7 is even, [/2] = [(m + 1)/2] and n — 2 [n/2] = 0 so that (2.6) holds for 
n+ 1. If nis odd, n — 2 [m/2] =1 and [(m + 1)/2] = [n/2] +10 that the last 
term becomes 
Crna] Frti- [(m +1)/ aif nies 2[(n+1)/: 2] glin+1)/2] ! 

which is the last term of (2.6) for 7+ 1. We note that the coefficient of cj 
remains unchanged and hence always has the value unity. 

Inserting (2.6) into (2.5), we have 
on Uy, 2) o" UR, 2) 


orn or” 


27 


a4] ai [c"(e) ao cosn fb Fe-) (2 r — 2 a cosd)"-*J 2) dh do] 


or 


n! u"(O, z) = n! u5(0, 2) 
1 22a 


+ _ ds 4 C; giao} hee c"(o) 0” -3+1 cos"-2J 6 cos nb F"-*J dd do . 
It 


But 
cos*d = aa cosa d + a, cos(a — 2)h+-:: 


so that the only term which will survive integration with respect to ¢ is the 
term 7 = 0. We have then 
1 
— 4)" 
na! wX(0, 2) = n! wh(0, 2) + <I" fore) ott Endo, n #0 
é 


and 
1 


U%0, z) = U8, 2) + 5 “i £%0) 0 F° do. 


0 


i tg 2} 4 


——- = 


But F is a function of 6 + 2 and so the result (I) follows. 
We remark that the apparent difficulty in (Z) for z = 0 is only illusory. 


For by lemma 1 
Qn etku 


Cr(o) = o" er(@) sothat o"(9) e"" Seam 


remains bounded for z = 0. This remark is important for our later results. 


ae eee KS eT 
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3. Axially Symmetric Solutions 


For »=0 and k=0 equation (2.4) becomes the equation for axially 
symmetric potential functions. It is well known that such functions are uniquely 
determined by their values on the axis of symmetry, provided those functions 
are analytic. In fact in terms of those functions the solution can be written 
down explicitly in the form, 


BANATERA EDF. (3.1) 


Vr — ¢ 0? 


U(r, 2) = 


These facts can be extended to the general equation (2.4),, as discussed by 
Henrict [4]. Again a solution is determined by its values on 7 = 0, provided 
these functions are analytic. Further, the solution can again be written down 
in the form, 
eae Dn A) Nee a : 2 2\n/2—1/4 v2 pd 
Bee 2) (2) ni ee naa (h 7? — 64) (3.2) 


= fh 


x w (0, 2+ 70) do- 


Equation (3.2) shows that actually a complete knowledge of w(0, o) is not 
necessary, but rather that only w(0, z+ 720) + w (0, z—70) need be known. 

We obtain our second fundamental relation by setting z = 0 in (3.2), thus 
obtaining 


w(r, 0) = 


P(m+1) (k\-n+2 1 ae 
Tug (2) wef — eV naa (bY — e) (1) 


0 


x [w(0, 49) + 0, -1 e)] dg. 


Formula (II)is aVolterra integral equation for the quantity w(0,70) + w(0,—i@). 
Its exact solution is not of concern to us as we shall encounter only the case 
in which the left hand side is zero, and then w(0,7@) + w(0, — i) will be 
identically zero. 


4. Solution for Small Wave Number 


In this section we shall derive integral equations of second kind analogous 
to those of Jones [3]. These are particularly well suited to approximations for 
small k. We first derive the equations formally and then show that they do 
indeed yield solutions of the problems. We content ourselves with the deter- 
mination of the cosine terms, that is, the functions U*(r, z), the extension to 
sine terms being clear. 


=_T AN SN eee 
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We begin with the formula (I) and write 22—=— «. Upon integrating (I) 
m times with respect to «, we obtain 


nt | -- [ [U-(0, iz) — Ux(0, iY z)] ar ++ dey 
i 6 (4.1) 
PAPE AT ek cd gee) | (8) —— 


where the A’s are constants and K,(8) = 0" c"(0), 0? = 8. This expression was 
derived under the assumption that « was real and negative. We observe, 
however, that all of the quantities which occur in (4.1) are analytic functions 
of the complex variable «, as long as « is not real and positive. If we choose 
the branch of VB — a which is positive for arg a = 0, % < f, the integral on the 
right side of (4.1) is analytic in the complex «-plane cut along the real axis 
from 0 to 1. Thus w*(0, i Vx) will suffer a discontinuity on approach to the real 
axis, 0 < « < 1 from above or below. Indeed from (4.1) we have for arg « = 0, 


n! = | [w(0. tyz) + ur(0, i¥r-) — ut(0, it) — un(0, i/2-)] de +++ de, 


0 


n—-1 
+27 A, a" = T(x; K,) . 
m=0 


Here w”(0, i Vz-) denotes the value which u"(0, Vz) assumes, when arg t = 0, 
after a clockwise circuit of the origin and 7"(«; g) denotes the operator, 


T"(a; g) 


“fa ) See ape | lh eth VB—a as! 


ae of jp ee | 
eS, |i = ip /ato ae ‘By 


Our next observation is based on formula (II). Since U = 0uoz=0,r< 1 
we have also w"(r, 0) = 0 on x < 1. Applying (II) to the functions u” we have 
accordingly, 


(4.3) 


u(0,ir)+u0,-—tn)=0, O<r<l, 


since these quantities will be solutions of a homogeneous Volterra equation, 
but this means that for arg t = 0 


w"(0, ijt) + u(0, iyt-) == (0), 


so that the first iterated integral on the left side of (4.2) is zero. The iterated 


(4.2) 
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integrals of wu” are known functions, so that except for the values of the con- 


stants /,, the left side of (4.2) may be-considered as known. For simplicity we 
write, 


at [oof ualo, iyt) + ut(0, iyt-) de ++ de, . = Ha) . 
0 0 


Now let us denote the two integrals on the right side of (4.3) by T7’(«; g) and 
T;'(a; g) respectively so that (4.2) becomes 


n—1 : 
H(a) + D' A, 0" = 6 T'(«; K,) + Te; K,) - (4.4) 
m=0 


Essentially the equation (4.4) for » = 0 was obtained by JonEs [3]. He made 
the observation that if one finds the inverse of the operator Tj'(«; M,) and 
applies it to (4.4) the product (7;’)-+ 77’ is an operator which is simple. It is 
not essential that one works out the exact details of this process but for com- 
pleteness we outline the details. If we write 6 = 1 — o anda = 1 — y, equation 
(4.4) becomes 


H(L~7) + YA, (1 —7)™ 


1 o 
_ (=1)" FE Ke ee sinhk /o — y d cosk Vy —o 
= all — 0) ——_————__ do + | K, (1 — 0) ——— do}. 
2 i ( ) eo o A ( oy) ye S 


Now we recognize that the last integral may be inverted with the aid of the 
Laplace transform to give us 


(—1)” Kt i rie (eee sinhhk(/y + Vo) sinhk(//y — Vo) | 
Genoff ie? ee ee 


ee ee : not F 
7% Peal ¥ yo—y ay He — 7) + 3 Aglt — 7) |. 


Upon writing o = y?, » = x2, we obtain 


+1 


VK — 9") + ra fe K, (Uta) eee 


4 Cae 


y ——— 
Sete | Oni 
wt y i Vy? — 42 


n—-1 
% He (1 — x4) + D2 A, (1 = 2°) ahd 


m=0 


, 


" 
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The essential fact to observe about (4.5) is that it is a Fredholm integral 
equation of second kind with a regular kernel which is small when k is small. 
The method of successive approximations will certainly apply and yield the 
solution as a power series in & for k sufficiently small. We remark that in the 
actual mechanical process of computing the successive terms one need not 
make the transformation to (4.5). Since the operator 7, involves 
sinhh Va — 6 
ya—B 
it will produce terms of one higher degree in & than 7. Clearly an iteration 
process may be based directly on (4.4). That is, if we assume 


oo 
ae pS m 
K,= bee k 
m=0 


then it is easily seen that the K,,,, will be determined by recursion in the form, 


nm 


SAA eae Dhae 
/ /i— = d B a. Biker m—1 K,,0) 


a 


or inverting this simple integral equation®) 


de Wie 
Kam(%) aaa ‘ da. / PAIS syed er Fg) dp 0 
0 


In the developments so far there has been no apparent difference between 
the various values of n. The complication which arises is, of course, contained 
in the determination of the constants A,,. For this we use the remark made 
in connection with formula (I) that the quantities c”(o) must vanish like o” when 


o > 0. Recalling that M.(6) = puto ( /) 


this means we must have K,() vanishing as 8" when f -> 0, that is, 


gw KulB) = 0 SE Lyn, ST (4.6) 
Thus we attack (4.4) as follows. Define functions P"(«), Q””(«) as solutions of 

H*(g) = I*(a; P), a@= Tr(a; QO), m=0,1,...,.m—-1. (4.7) 
Then choose constants A», Ay, ... , A,-1 so that, 


qm 


dam 


n—1 
P(x) + DS” Ap O° (a) == 0 at ag 0. re 0, ee 1 8) 
m=0 


5) This is easily done with Laplace transforms. 


ZAMP XI/17 
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It would follow then that 
n—1 
K, (a y= Pug )+ ae QO" (a) 


would be a solution of (4.4) satisfying the condition (4.6). 

The equations (4.8) constitute a system of linear equations for the A,,’s 
which will have a solution unless the determinant should happen to be zero. 
If this were so then the homogeneous system corresponding to (4.8) would 
have a nontrivial solution A®,..., A°,. Following the general verification 
procedure which we present shortly, one would find that if 


as =o SA G22 
then 
U(r, 0, 2) = M,(r, 2; c2) cosn 0 


is a solution of equation (1.1), satisfying a radiation condition and vanishing 
on z= 0,7 < 1. The uniqueness theorem for the problem of section (1) then 
states that U(r, 0, z) =0 hence c%(9) is identically zero and finally all the 
constants A®, must be zero. The contradiction shows that the determinant of 
system (4.8) cannot vanish. ; 

We turn now to the verification of the formal solution which we have 
derived. First we observe that the equations (4.7) do indeed have solutions at 
least for sufficiently small k. Now consider the function, 


n—1 
K,(a) = Pra) + DA, Orm(a) 
n=0 
with the constants A,, chosen so as to satisfy equations (4.8). We put 
c,(0) = o-” K, (0?) and consider the function, 
U(r, 6, 2) = [Uy (7, 2) + Mr, z;¢,)] cosn 0. 


U(r, 6, 2) will then be a solution of (1.1). Hence it will be an analytic function 
and lemma 1 applies to yield, 


Op 52) MAG C) = wee 


u"(r, z) being a solution of equation (2.4). But now uy, uv", and c, will be con- 
nected by the formula (I). Thus we find corresponding to (4.2) 


nl u"(0, iy) + wr(0, i fa) = na! wt(0, i Va) + un(o, iV") 


0 
a Our 


Vol. XI, 1960 On the Scattering of Waves by a Disk 259 


Note that the ~ derivatives of T(x; K,) may be formed even at « = 0 because 
we have constructed K,(«) to vanish like x" as « > 0. Now Pr(a) and Q"™(«) 
are solutions of the integral equations (4.7) with H”(x) representing the m-fold 
iterated integral of — [w%(0,i Va) + w3(0,iVo-)]. Thus differentiating (4.7) 
n-times and comparing with (4.9) we deduce, 


uO, i ya) + w(0, iPo-) = 0. (4.10) 


We have only to enter the result (4.10) in formula (II) to conclude that 
u'(r,0) = 0 for r<1. Thus U(r, 6,2) is the (unique) n-th Fourier cosine 
coefficient of the solution to the problem of section 1. 

We have thus exhibited a method for calculating approximately any finite 
number of terms in the expansion of the function U(r, 0, z) of section 1. The 
degree to which an actual solution is approximated will depend of course on 
the ease with which the incident field, U(r, 6, z), may be approximated by a 
finite number of terms of its Fourier series. 

It would perhaps be well to insert here a few remarks concerning the small k 
approximation. Let us suppose that the function H"(«) has the form 


co 
H*(a), = ge ee 
m=0 
Then we write 
fe 
Tia Se 
m=0 
What the iteration procedure for small & yields then are functions Py(«), Qi" (a) 


+N 
such that 
Ht — Tr(a; Pt) = O(N); a — Tr(a; OR") = O(R*}). (4.11) 


We choose constants A so that equations (4.8) are satisfied for the approxi- 
mating functions Py and Q%". Then set 


ch (0) = 9" Pyle?) + S' Ay, On" (4) 
m=0 
and consider 


UN (r, 0, 2) = (U(r, 2) + M,(r, 2; cf)] cosn 6. 
Again UY will be a solution of (1.1) so we can write, 
Un(r, 2) + M,(r, 2; of) = 17" unl, 2) - 
Then following our previous procedure we deduce that 


ut(0, iV) + ui(0, i Var) = 0(kY*) 
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which again by formula (II) yields 
U(r, 6,0) = 0(RN*") for r<l. 


Thus solution of equation (4.7) to terms in k* yields the Fourier coefficients 
of the function U(r, 6, 2) of section (1) to the same degree of approximation. 


5. Diffraction by a Hard Disk 


The problem considered by JoNEs [3] was the diffraction of a sound wave 
by a hard disk, that is, the condition (1.2) is now replaced by the vanishing 
of the z-derivative on the disk. In this section we indicate briefly how our 
methods need to be modified to handle that case. We shall maintain the same 
notation as in the previous sections, with U(r, 0, z) to satisfy 


(7,050) = 0 StOrae alas (5.1) 


instead of (1.2). Replacing (1.3) we have 


1 22 


Ulr, 0, 2) = Uylr, 0, 2) + z= e | [ted Re ydbde ae 


where /(g, ¢) this time represents the discontinuity of U itself on the disk. 
If we require that the solution be continuous it follows that 


f(o,¢) +0 as o->1. (523) 


The functions U, U, and f can be decomposed into Fourier coefficients as 
we have done in section 3. In place of (5.2) we get by an integration over z 


Zz 


2 1 2x ik : 
[00.6.9 Hr, 6,nare= | [Hedy eddae, 
Oe) 


0 
where c is an arbitrary constant. The results of section 3 now carry over with 


minor modification. For example, when r = 0 


1 2% 


Jivo.0.9—u, WOO Naren se | [fed edb de, 


where yu = 9? + 22. If we write 
2% 


[ U0, 6,4) d0= wy), / U,(0, 0, t) dO = u0(t) and [ile d) db = cy. 
2 0 


v 
¢ 


il 
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This last equation becomes 


2 


1 
[we — weolate= > [ele 
0 


0 


etku 
“ag. (5.4) 


We now continue equation (5.4) into the domain of complex variables and 
we find that if. z>7zy (0 < y < 1) we get 


ty ¥ ——- 1 —— 

“ 1 “ —kVy*—o2 x, ikVor—y 
| (w(t) — uf()] dé +e= 5 | / co(Q) @ === dp + / (leit 7 
0 0 es y 


while if z > —zy (0< y< 1) we get 


Y yh Estas hVe®aot : ikVor—y 
| (w(t) — up(Z) dé] + o = > -3| / — ¢4(@) ae do + tr (0) 9 vl. 
a ‘ ee j 


Upon adding these two equations and making some trivial changes of variables, 
we get 
> 


[cet t) — u°(—i t) — u9(i t) + uQ(—id]idi +2 


d : k Vy? a kVo—y 2 © Z 

° <sinhk /y?—o cosk Vg? — y® 

— = o)d do}. 
, ee ee) joom “wele | 


Now the boundary condition (5.1) implies, according to the representation (3.2), 
that w(t t) — w°(— 7?) is a constant c,. Hence equation (5.4) assumes the form 


y 
cy +2e + | (ue(-9 — u9(i t)] i dé 
0 


a ee : 
° 4sinhk /y? — cos k — 
/ 1sin Ms 0? 0 C(O ) do Jas — zi col 2a ; 


The constants c and c, are evaluated from the requirements that co(o) > 0, 
0 > 1 and ¢,(0) is finite at @ = 0. The equations involving the higher angular 
harmonics may be derived in much the same manner and indeed involve only 
further integration with respect to the variable z* = (-- y*) but we shall not 
pursue this matter further. 


6. Incident Plane Wave 


We include here some numerical results for the case of an incident plane 


wave for Dirichlet boundary conditions on the disk. There is clearly no loss 
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of generality in assuming that the plane wave is incident upon the disk from 
the side z < 0 and that its propagation normal lies in the x — z plane. Let y 
denote the angle which this normal makes with the positive z axis. The incident 
field in rectangular coordinates is then 


eth (xsin y +zcos y) 


and in cylindrical coordinates it has the form e**("sinvces?+#0087), Hence 
1 2% 


Ulr, 6, 2) 22 gik(rsinycos@+zcosy) 1 aa 4 0, 8 i do’ do. 


On the axis 7 = 0, this representation assumes the form 


1°27 


tkzcosy 1 ae fe es r 
U(0, 8, 2) = eter + Tf I flo, 0) 9.6 do. 
(ent) 


Or, upon integrating over 6 and 0’, we get 


1 


u9(0, 2) = eb PRcosy z= [%) s 


0 


iku 


odo. 


If we continue z into the domain of complex variables and use the fact that 
uX(0,7y) + u(—ty)=0, Om yal, 
we get 


c —— 
cosh (k y cosy) ale 2 . as yy 
; ye@—y 


t c%(t) dt 
(6.1) 


2 72 
if a Fe) aha Of 


ye — 72 


Upon rewriting equation (6.1) in the form which we described in section 4, we 
obtain 


1 
def sinhk (B +0) . sinhk (B — o) 
+= / 6B) | Ea uae “| ap 


5 (6.2) 
Oe Wk mer h (k Vo® — 
ar aa rd ‘| g do cosy (y/1 — 9 cosy) yee e 

0 


Here Go) = V1 — y2 cy) and o=V1—y%. Upon applying the Neumann 


age a “eo 
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NO 
j=) 
Oo 


iteration procedure we find for small k 


+ { 20k sin? 0 4 
Ores ee Wp oe Pre eeu 
c%(0) Af gece | - --k IS "4. @2 cos?y “al 
__ = pg [cosy i ae 8 
ah |“ 3a on ee y 
+ ae [SE? (Se = 40% = 1) 2s (lai OF a au) 
ee 1 jy 2 0 + 08) sels (6 — 2 0? + 2 cosy) 
24 a 0 T OQ Bee 0 4 COSs*yY 
+e ]} + 008). 


In the limit &, = k (normal incidence) this result coincides with the one which 
Bouwkanmp [1] obtained. 
For the first harmonic, we obtain the equation 


De eke inhk, y hk, 
—4iksiny [y ee Ff. a ») +A | 
1 ; : : / (6.3) 
-| are Ve — y? we #2 cl (t) dt + 7 Ai sin aie —_ we #2 cl(t) dt, | 
— y — 
while for the second harmonic, we get 
1 
z fs Age 
/ Fe se WE” 8 2(t) dt + i 7 oe a “2 2(t) dé | 
TB ee aoa (6.4) 
orp a eas 3y oe fn 3 oe | 4 A2 AB ye | 
Zz 4 z 


Here the A’s are constants of integration. We observe that for all three cases 
which we consider here, that is, the zeroth, the first and the second harmonics, 
we are always solving the same type of integral equation, the only change 
being in the form of the non-homogeneous term. Equations (6.2) and (6.3) are 
solved in the same manner which we solved (6.2) save for the fact that we have 
to evaluate the constants of integration. These are evaluated by the require- 
ment that c1(t) and c2(t) are finite at ¢ = 0. We omit the details and merely give 
_ the results. We have 


(o) = — StF SEY [1 + % (2? —1) cost y — (1 + orl} + 0184), 
ayi—¢ oe 
2 
(0) = SEE? + of). 


1 — @? 
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The terms are calculated to different orders of k so that it may be possible to 
calculate the scattering cross section to terms correct to the order #* if need be. 

The work on this paper was sponsored by the Office of Ordnance Research, 
US Army and the US Air Force Office of Scientific Research. 
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Zusammenfassung 


Es wird gezeigt, dass das Problem, die Beugung einer ebenen Welle an einer 
kreisférmigen Offnung oder Scheibe zu bestimmen, auf die Lésung von regularen 
Iredholmschen Integralgleichungen zweiter Art zuriickgefiihrt werden kann. Die 
Lésungen dieser Integralgleichungen liefern uns fiir die Scheibe im Falle der Neu- 


— 


mannschen Bedingung die radiale Variation der Unstetigkeiten der Wellenfunktion — 


und im Dirichletschen Falle die Unstetigkeiten ihrer normalen Ableitung. Ist das 
Produkt von Offnungsradius und Wellenzahl klein, so kénnen die Integralglei- 
chungen gelést werden. Fiir die Ableitung der Integralgleichungen verwenden wir 


einerseits die Poissonsche Darstellung fiir die Wellenfunktion und andererseits die — 


Fortsetzung der Helmholtzschen Darstellung in die komplexe Ebene. 


(Received: July 25, 1959.) 
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Solution of the Poisson Equation Related to a Boundary-Value 
Problem for the Sphere 


By Evan E. Jones, Nottingham, England!) 


Introduction 


There are three boundary-value problems associated with singly-connected 
regions fundamental to the theory of potential ([1], p. 213)2). The interior- 
type problems associated with the region V bounded externally by a closed 
surface S having continuous curvature are (i) the Dirichlet problem, in which 
a function, harmonic in V, is required, having prescribed continuous boundary 
values on S, (ii) the Neumann problem, in which a function, harmonic in V, is 
required, having prescribed continuous boundary values for its normal deriva- 
tive on S, (iii) the problem usually referred to as ‘the heat conduction problem’, 
where a function, harmonic in JV, is required, having continuous boundary 
values for itself and its normal derivative on S, which satisfy a given linear 
relation. In conjunction with these interior-type problems are the exterior- 
type problems where functions, harmonic in the infinite region external to S, 
are required, satisfying similar boundary conditions, with prescribed forms at 
infinity. A more general type of boundary-value problem is the exterior-inte- 
rior type, usually referred to as ‘the electrostatic dielectric problem’ ([2], p. 194), 
in which harmonic functions are required in regions external and internal to a 
surface S, satisfying two boundary conditions involving the continuity of the 
functions and a prescribed discontinuity in the normal derivatives on-S. 

A solution to the last mentioned problem has been determined by YEH et 
al. [3], when S is the surface of a sphere, in terms of harmonic functions having 
prescribed singularities in the regions external and internal to S. The solution 
is stated in the form of a sphere theorem, which is the culmination of less 
general theorems by various authors. The solutions of the first three boundary- 
_value problems for the sphere appear as particular cases of the general result. 

It is proposed in this paper to present a solution to a similar type of problem 
with controlling differential equations of the Poisson type. The range of phys- 
ical problems which can be solved is thus considerably extended. 


1) University of Nottingham, Dept. of Mathematics. 
2) Numbers in brackets refer to References, page 272. 
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The Solution of the Problem 


It is tequired to determine the functions @,(7, 9, @), y;(r, 0, w), having 
prescribed singularities, and satisfying the Poisson differential equations 


V? pe = flr, 8, @) , 
V29,=1: (7, 8,0), (2) 


in the regions V,, V;, respectively external and internal to S, the spherical 
surface of radius a, with centre at the pole of spherical polar coordinates 
(r, 0, w). The functions ¢,(7, 9, w), p;(7, 9, w) and their derivatives 0p,/0r, 0@;/07, 


normal to S, are required to satisfy the boundary conditions 


Pe = YP; ’ (3) 


OPe Og; 
cy Fe + B,D, = Hi; a +B 9: +Y.- (4) 


In (1) and (2) the functions f,(7, 0, w), /,(r, 0, w) are together with their first 


partial derivatives assumed to be continuous in the regions V,, V; respectively. — 
Also, in (4), the coefficients «,, 6,, %;, 8; are given real constants, and w(0, w) is — 


a prescribed continuous function on S. 

The solutions are developed in the forms y, = y,, + D,, 9; = Y,; + D;. Here 
Pe”, 9, w) {or y,;(7, 8, w)} is a particular integral of (1) {or (2)}, which together 
with its first partial derivatives are continuous in V,+ S {or V; + S} except 
at a finite number of singularities prescribed in V, {or V;}, and @,(r, 0, w) 
{or ®,(r, 0, w)} is a harmonic function, being single-valued and continuous in 
V,+ S {or V; +S}. The particular integrals m,, and g,; are either prescribed 
functions, or in general can be determined as the sum of two types of functions. 
Thus 


Moe al P. an OM Poi = Pr; ¢e Pri, where px.(7, 0, @). per, 6, ) 


are Newtonian potentials due to volume distributions of density —/,/4z, 
—},/4 in V,, V; respectively ([1], p. 150) of the form 


and Poel’, 0, w) is a harmonic function having the prescribed singularities in 
V, (including infinity), and g**(r,0,@) is a harmonic function having the 
prescribed singularities in V,. It is noticed that the functions Poe Po; and the 
normal derivatives d—,,/0r, Op, ;/Or are continuous on to S. 


The functions @%,(r, 0, o), p*,(r, 0,@) are readily determined in certain 


special cases. Thus if f,(r, 0, ) and f,(r, 0, w) are functions of one space coor- 


ee mr 
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dinate only the differential equations (1) and (2) can be integrated directly to 


yield particular integrals, e.g. when f, and /; are functions of the radial coor- 
dinate 7 only, then 


w= > io rar — ie rar, 


and 
ir . 
P= — > [hl Par — [fy rar. 
0 r 


On substituting for ¢,,(r, 0, w) and g, ;(r, 0, w), and their normal derivatives 
on S into the boundary conditions of (3) and (4), it is possible to write 


(D, — ,), = g(6, w) , (5) 
(a, Se + BP, — a, 2 — B,,) = h(6, 0), (6 
where by definition 
(Poi — Poe)s = 8(9, @) , (7) 
(a; Fe + Bi gys +p — te Pee — Bg.) = (0, 0) 6) 


The functions g(8, w), 4(@,@) are continuous on S, hence by use of a known 
solution of the Dirichlet problem it is possible to form two unique functions 
g,(7, 9, w), h;(r, 0, w), which are continuous in V,; + S and harmonic V;, having 
(0, w), h(6, w) respectively as boundary values on S. The functions are deter- 
mined in the form of Poisson integrals ([1], p. 240), such that 


ae) 0’, w’) 
g:(7, 6, w) so a = + Bt a : aS , (9) 
$ 
and ; : #3 
h,(r, 0,0) = 4 —* ib EGS, (10) 


where R is the distance from the point (7, 8, w) in V; to the point (a, 6’, w’) on S. 
Y The harmonic functions ®,(r, 0, w) and ®,(r, 0, w) occurring in the complete 
solutions for y, and y; respectively are determined by the equations 


(a, + a) (7) = 2% 9 (S) — 2 (FS), (11) 


td Y id a 


(a; + a.) P(r) = —&%, g(r) — U(r) , (12) 
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where 
1 
I(r) = [{ahilrt) +» g(r} eat, (13) 
F 
in which 
he _ ee Be) J (14) 
Wess on as As : (15) 


and for simplicity of expression g,(7) =g;,(7,9,@), etc. If g,(r,@,m) and 
h,(r, 9, o) are harmonic functions defined in V, with boundary values g(0, ), 
h(8, wm) on S, then by the Kelvin inversion transformation ([1], p. 231), it is 
seen that 

a? 
Y 


g(r, 8, 0) =~ 9(<, 6,0), (16) 


oe a 
h,(r, 0,0) == = n(—, 6, o). (17) 
In general if g;(r) ~ 0(1) for ry > O, then g,(r) ~ O(1/r) for > oc. Similarly for 


the function pair h;(r) and h,(rv). Alternative forms for the solutions of (11) and 
(12) exist as 


(1; + a4) P(r) = a4 g(r) = Lr), (18) 
(a; + a) ®,() = — 8% g(@) — 2 a(S), (19) 

where 
I(r) = Je hr) ty gla} en dt, (20) 


and thus J,(r) = a/r I;(a?/r) by comparison of (13) and (20). 
These results hold in the given forms in general when mu > 0, and in partic- 
ular for w = 0, providing 


kc h(0, 0) +9 g(0,a)} dS =0. (21) 


On integration by parts it is possible to extend the range of existence of the 
solutions to negative values of “, with the exception of the negative integers 
= —1, —2,.... Thus, for example, considering the integral 


a 


iG t) "dt 


0 
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of (13) defined for u > 0, it is possible to integrate by parts leading to the 
expression 
1 


1 an a 
(el ae "ly #) pee 
: (7) 7 é h,(r t) t* dt , 


and consider this as the continuation of the original integral to the range 
—1< y4< 0. Repeated integration by parts will further extend the range of 
application. 


The Verification of the Solution 


In the first instance it is necessary to prove that @,(r) is a harmonic function 
without singularities in V;, and that @,(r) is a harmonic function without 
singularities in V,, being regular at infinity. 

The Poisson integral g,(r) of (9) is continuous in V; + S and harmonic in 
V;, having the development g,(r) = g; (0) + O(r) for ry <b < a. Thus (1/r) g,(a?/r) 
is continuous in V, + S and harmonic in V,, decreasing to zero at infinity to the 
order O(r~4). Similarly for /;(r) and (1/7) h; (a?/r). 

The integrand of J;(r) of (13) is harmonic in V; and continuous in V, + S 
for 0 <¢<1. The integral exists and is uniformly convergent for 7 < a if 
fu = 0, with the restriction that ah;(0) + g,(0) =0 when uw =0, which is 
consistent with (21). By application of a modified form of Harnack’s first theo- 
rem ({1], p. 248) it follows that J,(r) is harmonic in V;, continuous in V; + S. 
Similarly J,(a?/r) is harmonic in V,, continuous in V,+ S, and decreases to 
zero at infinity to the order O(r-). 

From (11) and (12) it is thus seen that @,(r) and D,(r) have the required 
harmonicity properties in V, and V, respectively. 

Secondly it is required to prove that @,(r) and @,(r), and their normal deri- 
vatives on S, satisfy the boundary conditions of (5) and (6). In this connection 
it is required firstly to establish that these functions and derivatives are 
continuous on to the boundary S attaining finite limiting values there. 

Since g;(7) can be written as the sum of the potentials of a simple layer and 
_ double layer distribution with continuous density proportional to g(6, «), then 

g,(r) is continuous on to S with a limiting values g(0, w) ([4], p. 245). Also the 

normal derivative g;(r) attains a limit on S as y >a in the region VV, if the 

limits of the normal derivatives of the above-mentioned potentials exist. A suf- 

ficient condition for this is that the function g(6, @) is not merely continuous, 
put has continuous first and second derivatives with respect to the space 

coordinates ([4], p. 151). Similarly for the function /,(7) and normal derivative 
Ayr). 

Lia I(r) is uniformly convergent in V,+ S for w = 0 it follows that J,(7) 
is continuous in this closed region, and in fact attains the limiting value J;(a)- 
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on S. Again as a consequence of Harnack’s first theorem it is possible to form 
I(r) by differentiating J,(7) with respect to 7 under the integral sign. On 


integration by parts it can readily be shown that 


I(r) = = {ahr +98} — 2 Le) 
for r <a, u = 0. Since all functions on the right hand side of this equation are 
continuous on to S, then so is J; (7). 

A similar argument shows that g;(a?/r), h,(a?/r), I;(a?/r), and their radial 
derivatives, are continuous on to S as 7 > a@ in the region V,,. 

On substitution for ®,(7), ®;(r), and their radial derivatives, from (11) and 
(12) into (5) and (6), and allowing r > a, it can be verified that the boundary 
conditions of (5) and (6) are satisfied. ; 

It has thus been shown that the expressions (11) and (12), or (18) and (19), 
lead to solutions of the boundary-value problem. 


Arbitrariness in the Choice of Particular Integrals 


The particular integrals @,,(r) and g, ;(v) of the differentia] equations (1) and 
(2) consist partly of harmonic functions y**(r), p**(r) having prescribed singu- — 
larities in V, and V; respectively. It is required to consider the case with 7 
prt (r), pr*(r) replaced by ga*(r) + 2,(r), pt*(7) + Q;(r), where 2,(7), 2,(r) are 
harmonic functions without singularities in V,, V; respectively. On forming the 


altered values of g(9, w) and h(6, w) from (7) and (8), it is possible to form the 
functions 


Y 


G.(r) = gi(r) + 2) —2 Q(Z), 


Hyp) = hip) + a, 5 O(F) — Fe" O(F) + 2". 20). + 8,247), 
h, 


Y a or 
where G;,(r), H,(r) are the altered values of g(r), h;(7) respectively. These funa 
tions are harmonic without singularities in V, having the required prescribed 
boundary values on S. They are thus uniquely determined. 
If the altered values of ®,(7), ®,(r) are O*(r) and OF (r) respectively, then 
from (11) and (12) it can be shown that 


Since the altered values of ,(r) and g,(r) are ,,(r) + 2,(r) + @*(r) and 
Poi(”) + 9;(r) + ®¥(r) respectively, then the addition of a function harmonic 
without singularities in V, {or V;} to the particular integral @,,(r) {or @,,(7)} 
will not affect the value of the solution 9,(7) {or y;(r)}. The choice of particular 
integral is arbitrary to this extent — a suitable choice may however simplify the 
‘analysis leading to the solution of the problem. 
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Solution in Terms of Prescribed Singularities 


The prescribed singularities are included entirely in the parts p**(r) and 
Poi (”) of the particular integrals ~, ,(r) and gy, ;(r). Using the symbol (*) for all 
functions corresponding to the parts y*,(r) and p*;,(7) of the particular integrals 
it is possible to form 


a iq? 
gilr) = et (r) + * se tala 
ar Opety ro) 2 F 2 
hi(r) = hi (r) — 2”. SP880) _ gente) — a, 2 gte(S) 4 2% gax(*). 


The functions g;(r), ;(r) are harmonic without singularities in V,, and on S 
are consistent with the boundary conditions of (7) and (8) respectively. They 
are thus uniquely determined, and on substitution into (11) and (12) yield the 
results 


P(r) = DF (r) + yl), Di(r) = OF (7) + vil) , 


where y,(r) and w,(7) are functions identical (but for a change in notation) with 
the perturbation potentials of YEH et al. [3]. 


Some Particular Cases 


The examples listed in this section refer to problems controlled by diffe- 
rential equations of the Poisson type. For the internal-type problems it is 
assumed that (3) and »,(r) are redundant, and for external-type problems (3) 
and ¢,(v) are redundant. 

(a, = 6, =a;=0, B; = 1] Internal Dirichlet problem with ®,(r) = —h,(r), 
and (;)s = — (6, ). | 

faps= 2; = «, = 0,8, = 1) External Dirichlet problem, with (7) =H, (7), 
and (g,)s = (0, ©). | 

The solutions of the Dirichlet problem are uniquely determined ([4], 
pp. 183, 186). They have application to the electrostatic field theory for coa- 
ductors in the presence of volume charge distributions, and heat conduction 
with distributed heat generation and prescribed surface temperature on S. 

[a, = B, = B; = 0, «; = 1] Internal Neumann problem, with uw = » = 0. The 

tunction @,(r) is given by (12) with the restriction 


J h(9, 0) dS =0, and (28) = —y(0, w) . 
S 


(x; = B, = 8; = 0,«, = 1] External Neumann problem, with » = 0, oe =I. 
_ The function @,(r) is given by (11), with (0g,/0r); = y(8, w). Pe 
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The solution of the internal Neumann problem is uniquely determined to 
within an additive arbitrary constant, and the external problem has an unique 
solution ([4], pp. 183, 186). Application is possible to heat conduction problems 
involving distributed heat generation with prescribed surface temperature 
flux on S. 

lop = 8B, = 0, 0, =1, 27— 41 Internal third boundary-value problem, with 
vy = 0, 4 = ah. The boundary condition is 0y,/0r + hp; = — y(8, w), and ® (7) 
is given by (12). 

[fo 6,= 0, a,—1,.8, h] External third boundary-value problem, with 
y=0, 4=1+ah. The boundary condition is 0p,/0r—hy, = p(9,), and 
@,(r) is given by (11). 

The third boundary-value problems have unique solutions if h > 0. The 
problems apply to heat conduction with distributed heat generation under a 
surface condition with temperature determined by Newton’s Law of Cooling. 
The constant / is the radiation coefficient. 

[8; = B, = 0] These problems are of the internal-external types, and apply 


to dielectric potential problems with distributed volume charge, magnetic — 


problems with distributed permanent and induced magnetism, and heat con- 


duction problems involving two media of different conductivities in intimate — 


contact over the spherical surface S ((5], p. 17). Here v = a; wu, w = «,/(a; + &,), 
and since in practice w > 0 the uniqueness of the solutions is secured. If 
a“, =a, = 1, solutions are obtained for gravitational problems involving distri- 
buted mass. ; 
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Zusammenfassung 


___ Es werden Lésungen der Poissonschen Differentialgleichungen in den Bereichen 
innerhalb und ausserhalb einer Sphare, die gewisse Grenzbedingungen an der Flache 
der Sphare erfiillen, gesucht. Solche Lésungen werden schliesslich als Funktion 
partikularer Integrale der Differentialgleichungen ausgedriickt. Ihre Anwendung 
auf die theoretische Physik wird anschliessend behandelt. 
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Zur Thermodynamik teilweise wisstiger Hagelkorner 


Von RoLanp List, Weissfluhjoch-Davos!) 


1. Einleitung 


Erste Experimente im Hagelversuchskanal [1]?) deckten ein Phanomen auf, 
das scheinbar nie eingehender beobachtet und bei allen bekannten Hagel- oder 
Vereisungstheorien nicht in Rechnung gestellt wurde: das an bestimmte atmos- 
pharische Bedingungen gebundene Auftreten von Eiswassergemischen als An- 
lagerungen an bereits vorhandene Eispartikeln oder an kiinstliche Objekte (2). 

Dieser Effekt ist an einen ungeniigenden Warmeaustausch vom vereisenden 
Objekt an eine umgebende, kiihle Atmosphare (Temperatur unter 0°C) ge- 
bunden, der es nicht gestattet, durch Abtransport der Gefrierwarme alles an- 
fallende, unterkiihlte Wasser zu verfestigen. Die Folge davon ist, dass jeweils 
iiberschiissiges Wasser auftritt, welches unter Umstanden durch den Luftstrom 
mitgerissen wird oder sonstwie abtropft. Dies kann manchmal auch beobachtet 
werden, der Schluss jedoch, dass es sich dabei um alles nichtgefrierende 
Wasser handle (LuDLAM [3)), ist irrig. 

-Die Tatsache, dass sich bei gewissen Verhaltnissen auf einem vereisenden 
Objekt eine Wasserhaut bilden kann, die von innen her so weit gefriert, bis die 
Objekttemperatur 0°C erreicht, fiihrt zu einem ganz speziellen, im folgenden 
naher erlauterten Gefriermechanismus. Die Temperatur an der innern Grenz- 
flache der Wasserschicht ist im stationaéren Fall aus thermodynamischen 
Griinden auf 0°C (siehe Figur 1). An der Grenzflache Wasserhaut — umgebende 
Luft wird das Wasser jedoch unterkiihlt. Das Gefrieren kann von aussen her 
nicht stattfinden, da normalerweise entsprechende Eisbildungskerne fehlen. 
Zusatzlich zum Temperaturgradienten, beeinflusst durch die Luftstromung, 
treten nun wirbelartige Zirkulationen in der Wasserhaut selber auf — vielleicht 
kann man sogar von Konvektionszellen sprechen — die das unterktihlte Wasser 
értlich zur Eis-Wasser-Grenzflache transportieren. Dies erméglicht ein schlag- 
artiges Wachstum von Eislamellen von innen nach aussen. In diesem Momente 
ware auch ein Gefrieren der 4ussern Begrenzung der Wasserhaut moglich — es 
kommt jedoch nicht dazu, wenn sich der Nachschub unterkiihlter Tropfen als 
_ Anlagerung an die Wasserhaut stetig fortsetzt. So kénnen wir uns vorstellen, 
dass ein Wachstum eines Eiswasseransatzes moéglich ist, indem sich die Was- 
serhaut jeweils vom Zentrum unseres Objektes weg fortpflanzt und dabei vom 


1) Forschungsstelle der Eidgendssischen Kommission zum Studium der Hagelbildung und der 


Hagelabwehr am Eidgendéssischen Institut fiir Schnee- und Lawinenforschung. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 305. 
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bestehenden Eis oder Eisgeriist her durch Eislamellen oder Eiswande stabili- 
siert wird. Die Dicke derselben hangt ab vom Betrag der Warmeabfuhr nach 
aussen sowie den durch die Wasserhaut eingefangenen, unterkiihlten Tropfen, 
das heisst vom totalen Warmeaustausch Objekt - Umgebung, der den maxi- 


Konvektions - 
zellen 


wachsendes 
Eisgerust 


unterkuhite 


Wasserhaut 
Wassertropfen 


ta: 


Figur 1 
pes der Entstehung eines Eiswassergemenges als Ansatz an eine Eiskugel. tg Temperatur 
er Eiskugel, tx Temperatur der Grenzflache Wasserhaut/Eiskugel, ty, Temperatur an der 
Grenzflache der Wasserhaut gegen die Luftstrémung. 


malen Betrag gefrierbaren Wassers bestimmt. Auf diese Art entsteht ein meist — 


zusammenhangendes Kavernensystem, das mit Wasser aufgefiillt ist. Schleu- 
dern wir die fliissige Phase in einer Zentrifuge bei 0°C heraus, so bleibt nur das 
Eisgeriist tibrig. Dessen Struktur ist abhangig von den Entstehungsbedin- 
gungen und ist fiir einen bestimmten Fall in Figur 2 dargestellt. Die hellen 
Zonen bedeuten dabei durchsichtiges Eis, wahrenddem die eingebauten Luft- 
einschliisse durch ihre Randzonen dunkle Schattierungen ergeben, die einander 
tiberschneiden und tiberlagern kénnen. . 

Die Entdeckung der Eiswassergemenge iiberholt somit alle Annahmen 
und Berechnungen iiber maximalste Vereisungen (SCHUMANN [4] und LupLAM 


a 
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Figur 2 
Eisstruktur eines Eiswasseransatzes an eine Testkugel (d = 2. cm) nach dem Herausschleudern des 
die Zwischenraume des Eisgeriistes ausfiillenden Wassers. Entstehungsbedingungen: Lufttemperatur 
tr = — 30°C, Luftgeschwindigkeit v = 24 m/s, freier Wassergehalt w, = 8,1 g/cm’, resultierender 
Wasseranteil W,, = 19% des totalen Ansatzes (M 1: 0,13). 


[3], [5]), die auf der Basis einzig méglich erachteter «reiner» Vereisungen ohne 
Einbau von fliissigem Wasser beruhten. 

Das Ziel dieser Arbeit ist es nun, fiir das Zustandekommen von E1swasser- 
gemengen verantwortliche atmospharische V erhdlinisse auszumessen und sie in Zu- 
sammenhang zu bringen nut dem resultierenden Wasseranteil des totalen Ansatzes. 
Es sollen auch Formeln entwickelt werden, die diese Relationen beschreiben. 
Die Resultate werden zeigen, welche enorme Bedeutung dieser speziellen, neu 
entdeckten Form atmospharischen Eiswachstums bei der Bildung von Hydro- 
meteoren und auch bei der Erfassung des realen Geschehens bei Flugzeugver- 


eisungen zukommt. 


2. Theoretische Betrachtung des Warme- und Stoffaustausches 
bei Eiskugeln 


Zum besseren Verstandnis der thermodynamischen Grundlagen, die zu Eis- 
wassergemengen fiihren, seien die in Betracht zu ziehenden Vorgange erlautert 
und versucht, sie durch entsprechende Forineln zu beschreiben. 
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Gegeben sei eine Atmosphare mit der Temperatur ¢, (f, < 0°C!), dem 
Druck p, dem Wasserdampfgehalt wp, dem freien Wassergehalt in Form 
unterkiihlter Tropfen w, sowie deren Gréssen und Gréssenverteilung. Die 
Wassertropfen sollen dabei die gleiche Temperatur haben wie die Umgebungs- 
luft! In dieser Atmosphare werde zur Zeit t = 0 eine Eiskugel eingefithrt mit 
dem Durchmesser d, der Dichte 9x und der Relativgeschwindigkeit v gegentiber 
der Umgebung. Die Testkugeltemperatur ¢, sei 0°C. Diese letzte Annahme sei 
deshalb eingefiihrt, um von allem Anfang an eine Aufnahme von Gefrierwarme 
durch die Ursprungspartikeln zu verhindern. 

Mit der Zeit vergréssert sich unsere primare Eiskugel, der Ansatz kann unter 
Vernachlassigung einer méglichen Verdampfung durch folgende Gleichung 
charakterisiert werden: 


dm= E F,vw,dt —dm’'. (1) 

Dabei bedeuten 

dm Menge des an die Testkugel angelagerten Wassers [kg] 

E Einfangswahrscheinlichkeit fiir die an die Testkugel angeschwebten 
Tropfen [—] 

Fx Querschnittsflache der Testkugel senkrecht zur Bewegungsrichtung 
[m?*] 

v Relativgeschwindigkeit Kugel - Umgebung [m/h] 

Wy freier Wassergehalt der Atmosphare [kg/m*] 

at Zeitintervall [h] 


im Zeitintervall dr von der Testkugel wieder an die Umgebung abge- 
gebene Tropfenmenge [kg] 


Die Leistungssumme lautet fiir den gleichen Fall 


Qi = Okt + 2, te Ore (2) 
wobei 
OF; total weggefiihrte Leistung [kcal/h] 
Qe durch Konvektion und Warmeleitung weggefiithrte Leistung [kcal/h] 
Oy durch Verdampfung von H,O an der Oberflache des Testobjektes 
weggefiihrte Leistung [kcal/h] 
0, durch die angeschwebten, unterkiihlten Wassertropfen zugefiihrte, 
gegentiber der 0°C-warmen Testkugel iiberschiissige Kalteleistung 


{kcal/h] 


Diese total weggefithrte Leistung Q,,, kann unter der Voraussetzung, dass der 
Ansatz dm|dr gross genug sei, dazudienen, das angelagerte Wasser ganz oder 
teilweise zu gefrieren. Dabei ergibt sich die folgende Gleichheit: 


Ome = Qior (3) 


a 
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OnE Gefrierleistung der Menge dés pro Zeiteinheit angelagerten Wassers, 
das effektiv gefrieren kann [kcal/h] 


K6nnte alles angelagerte Wasser gefrieren, so wiirde dies die Leistung 
Qm, tot DenStigen, die sich ausdriicken lasst durch 


Qm,tor = (E Fy v wp — dm'|dt) S (4) 
mit 
S  Gefrierwarme von Wasser [kcal/kg]. 


Der Wasseranteil W,,[°%]| des Eiswasseransatzes betraigt damit 


Ar ae OmE als Qrot 
W,, = 100 (1 — ied, = 100 (1 — ee a (5) 
Es gilt nun noch, die weggefiihrten Leistungen Qx;, Qy und Q 
und deren Summe Q,,, in Gleichung (5) einzusetzen. 

Die vorteilhaftesten und vollstandigsten empirischen Gleichungen, die so- 
wohl den Warme- wie den Stoffaustausch von Kugeln beschreiben, diirften 
diejenigen von RANz und MarsHALt [6, 7] sein. Sie sind auch wesentlich all- 
gemeiner als die von KRAMERS [8], die jeweils in der Meteorologie verwendet 
werden (siehe Mason [9] und LupLaAm [3}), da sie den «Ventilations»-Koeffi- 
zienten genau charakterisieren und dadurch die Berechnungsmdglichkeiten 
fiir andere Drucke und Temperaturen geben. Die Gleichungen lauten (gemass 
Formeln 21 und 22 der zitierten Arbeit [6, 7]) fiir den in den nachfolgend 
beschriebenen Experimenten durchlaufenen Bereich der Reynolds-Zahl 
5-103 < Re < 5-104, unter Vernachlassigung der hier unbedeutenden addi- 
tiven Konstanten 


zu berechnen 


w 


Nu = 0,60 Pri? Ret? | (6) 
Nu’ = 0,60 Sct8 Rell? . (7) 
Dabei sind Nu Nusselt-Zahl 


Py Prandtl-Zahl 

Re Reynolds-Zahl 

Nu’ Analog-Nusselt-Zahl 
Sc Schmidt-Zahl 


In der Strémungs- und Warmelehre wird der Warmeiibergang durch Kon- 
_-vektion und Leitung mit Hilfe der Nusselt-Zahl Nw charakterisiert, wahrenddem 
- die Analog-Nusselt-Zahl Nw’ den Stoffiibergang kennzeichnet. Die Definitions- 


278 Rovanp List ZAMP 


gleichungen fiir die beiden Kennzahlen lauten 


Nu = shes, (8) 
Nui = ££. (9) 


Hierbei gelten 


Warmeiibergangszahl [kcal/m*h °] 
Kugeldurchmesser [m]} 

Warmeleitzahl [kcal/m h *] 
Stoffiibergangszahl [m/h] 

Diffusionszahl von Wasserdampf in Luft [m?/h] 


wrD~> AR 


Meist ist es jedoch vorteilhafter, den Stoffiibergang anstatt tiber den Kon- 
zentrationsunterschied Ac tiber den Dampfdruckunterschied des itibergehenden 
Stoffes zu definieren. 

Daraus resultiert die Gleichheit 


pAc=k,, Ap (10) 

mit 

Ac Konzentrationsunterschied des tibergehenden Stoffes langs des Trans- — 
ferweges [kg/m?] 

Ren Massenitibergangszahl {kg/m*h mbar] 

Ap Dampfdruckgefalle des iibergehendes Stoffes langs des Transferweges 


[mbar] 


Unter Beniitzung der Zustandsgleichung fiir ideale Gase kénnen wir die Kon- 


zentrationsdifferenz Ac auch wie folgt ausdriicken: 


Atc= LA : (11) 
yaa 
Dabei gelten 
R allgemeine Gaskonstante 
16 mittlere, absolute Temperatur im Transferweg des iibergehenden 
Stoffes [°K] 
M Molekulargewicht des iibergehenden Stoffes [kg] 
Gleichung (9) geht damit unter Beriicksichtigung von (10) tiber in 
, jer 
Nw’ = Mk (12) 


Oo Rowe. ~ 
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Ranz und Marsua tt [6, 7] gingen noch einen Schritt weiter und driickten das 


Produkt (R. T) wiederum nach der Zustandsgleichung fiir ideale Gase durch 
die Eigenschaften der Luft im Transferweg aus und erhielten so ihren endgiil- 
tigen Ausdruck fiir die Analog-Nusselt-Zahl Nw’: 


Nu’ = 2 ie (13) 
mit 
M Molekulargewicht der Luft [kg] 
p Mittelwert aus der Differenz Totaldruck minus Partialdruck des dif- 
fundierenden Stoffes im Transferweg [mbar] 
OL mittlere Dichte der Luft im Transferweg [kg/m*] 


(In der Definitionsgleichung von Ranz und MarsHALt fehlt die Grésse M im 
Nenner von (13), ihre Rechnungen sind jedoch korrekt durchgefiihrt.) 

Auf Grund der Gleichungen (6) und (8) lasst sich die Warmeiibergangszahl 
(durch Konvektion und Leitung) «,; wie folgt ausdriicken: 


ayy = 0,60 4 Pris Red? (14) 
Die durch Warmeleitung und Konvektion weggefiihrte Leistung wird damit 
Oxy = 0,60 Og 4 At Pris Ret? (15) 


At sei dabei die Differenz aus der Kugeltemperatur und derjenigen der Umge- 
bung, das heisst At = —t,. Wenn wir die Oberflache der Kugel Ox [m?] ersetzen 
iiber deren Durchmesser d, die Prandtl-Zahl Pr durch den im betrachteten 
Druckvariationsgebiet (1000 bis 200 mbar) praktisch konstanten Wert 0,71 
und die Reynolds-Zahl Re durch die iibliche Definitionsgleichung, so gelangt 
man zu folgendem Ausdruck: 


Oz yy 1,68 Ay (—t,) q3l2 yll2 | (16) 


Die kinematische Zahigkeit »[m?/h] kann dabei auch ausgedriickt werden 
durch y/o,, wobei 4 die dynamische Zahigkeit in [kg/m h] bedeutet. Letztere 
Grésse ist druckunabhangig bis zu freien Weglangen, die vergleichbar werden 
mit den vorkommenden geometrischen Dimensionen. Die kinematische Zahig- 
-keit » hangt somit vom Druck so ab wie der reziproke Wert der Luftdichte @,. 
_ Ebenso wie 7 ist auch die Warmeleitzahl 4 druckunabhangig in dem in der 
- Wolkenphysik zu beachtenden Druckvariationsgebiet. Es scheint nun vorteil- 
~haft, alle Gréssen bis auf den Partikeldurchmesser d und die Geschwindigkeit v 
_in eine einzige Druck-Temperatur-Funktion (P7),; Zusammenzufassen gemass 
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folgendem Ausdruck: 
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(PRee 1,684 Eee 


Diese Funktion ist in Figur 3 dargestellt. 


(PT) | ] i F 
[kcal/m? h””] Pig 
Rarer (P T)y 
6 : a | 
ee 
4 ' a 
“7 | eee —— eM; 
2 —_ ns == =e _ 
Dae Ne ae 
7 ES 
, as 
4 Ee 
IZ ZA f, 
-10 -20 -30 C 
Figur 3 
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Druck-Temperatur-Funktion (PT)x ,, (PT)y und (PT)xry in Abhangigkeit der Lufttemperatur 


tr, bei verschiedenen Luftdrucken p;. 


Aus Gleichung (16) wird damit endgiiltig fiir den Warmeiibergang durch — 


Konvektion und Leitung 


Qxr = (PT )g, PP vw. 


(18) | 


Wahrend der Stoffaustausch als Massendnderung nicht von Belang ist fiir die 
in Betracht gezogenen Verhiltnisse, spielt die dabei in Aktion tretende totale 
Verdampfungskalte L [kcal/kg] eine Rolle. Es wird daher vorerst aus den Aus- 
driicken (7) und (13) die Massentibergangszahl k,, berechnet. 


k,, = 0,60 


ho,M 
reo. 


Scll8 Retz, 


(19) 


: 


} 


Die bei diesem Stoffiibergang erzeugte totale Verdampfungskalte errechnet | 


sich nach dem Ausdruck 


oder 


Oy = km Ox Ap L 


07 = 0,60 aes L Og Ap Scil8 Rei, 


(20) 


(21) 


Unter L verstehen wir dabei die Verdampfungskalte, die aus dem Phaseniiber- 
gang fliissig-dampfférmig resultiert. Diese Annahme ist nicht ganz korrekt, da 


f 


| 
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nicht die ganze Oberflache einer durch Akkreszenz unterkiihlter Wassertropfen 
wachsenden Kugel nass zu sein braucht. Der damit verbundene Fehler, der in 
die massgebende Formel (26) eingeht, ist jedoch vernachlassigbar. 

Wenn wie vorhin bei (16) wiederum alle Gréssen bis auf die Faktoren Ge- 
‘Schwindigkeit und Teilchendurchmesser zusammengefasst werden, so entsteht 
eine 1.eue Druck-Temperatur-Funktion (P7),,, unter deren Beriicksichtigung 
sich ein. durch Stoffaustausch ausgelister Warmetransport ergibt von 


aes Q,, = (PT), a3 yi (22) 
mit ; 


we 
(PT),, = 599 sa Ap. (23) 


Gleichung (23) is’ ebenfalls in Figur 3 dargestellt. Im Gegensatz zu (PT) x; 
nimmt (P7), mit ©*igendem Druck ab, das heisst wahrend der Warmeiiber- 
gang durch Leitung und Konvektion mit abnehmendem Druck infolge der Re- 
duzierung der Molekiilaahl sinkt, wird ein Verdampfen von H,O-Molekiilen an 
der Partikeloberflache be. gleichen Vorgang begiinstigt und damit auch der zu- 
gehorige Warmeiibergang. Wie die Summierung der beiden Druck-Temperatur- 
Funktionen zeigt, ist diesy= Betrag im betrachteten Temperaturvariations- 
bereich bei Drucken zwisch»”i 1000 und 200 mbar praktisch gleich. Dass mit 
sinkendem Druck die Verdampfung den wesentlichen Beitrag zum totalen 
Warmeaustausch leistet, geht auch speziell aus Figur 4 hervor, wo das Verhalt- 


30 


(PTly _ Oy 
PT) y Oxy 
| in 
20 u 
Lnts, 
| 400 mbap 
788mbap 
| io00mbar 
=10 =) -30 ¢ 


Figur 4 
Verhaltnis der durch Verdampfen abgefiihrten Leistung Qy und der durch Konvektion und Warme- 
leitung abgefiihrten Leistung Qxz in Funktion der Temperatur tz, bei verschiedenen Luft- 
drucken py. 


“nis (PT);|(PT)x, gegen die Temperatur aufgetragen ist. Als Parameter 
wurde hierbei der Druck gewahlt. Diese Figur zeigt auch, dass der Anteil der 
durch Verdampfen produzierten Leistung Qy gegeniiber der durch Warme- 
‘leitung weggefiihrten Leistung Qx, mit sinkender Temperatur abnimmt. Ver- 
~dampfung spielt also bei steigender Temperatur eine gréssere Rolle. 
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Noch fehlt uns eine letzte Grdsse, um die Energiebilanzgleichung (2) explizit 
darzustellen, die Kalte Q,,, die durch das angeschwebte, unterkiihlte Wasser 
der Testkugel zugefiihrt wird: 


Oy = EF Oy (— ty) wy? (24) 
oder ¢ 


0, = 0,785 c, Bit) a, a i (ay 
c,, Spezifische Warme von Wasser [kcal/kg*].  _ 


w 


Die Einfithrung einer Temperaturfunktion bietet in diesem Falre keinen Vorteil. 


Gi "a \ ? 
—— ANG ry 
——-t=-10 °C 4 


Figur 5 


Beitrage zur Vernichtung der Gefrierwarme des angelagerten Wassers durch Konvektion und 
Warmeleitung Qx,, Verdampfen Qy und zugeftihrtes, unterkiihltes Wasser Q,, in Funktion des 
Wassergehaltes w,; v = 12 m/s, d = 2cm, ty = — 10°C bzw. — 20°C. 


. 

Um eine Vorstellung iiber das Verhaltnis der verschiedenen Beitrage zur 
Vernichtung der Gefrierwarme des angelagerten Wassers zu bekommen, seien 
in Figur 5 die Gleichungen (18), (22) und (25) dargestellt in Funktion des 
freien Wassergehaltes bei einem Kugeldurchmesser von 2 cm, einer Luftge- 
schwindigkeit von 12 m/s sowie den Lufttemperaturen — 10° und — 20°C. 
Figur 6 zeigt Verhaltnisse in Funktion der Luftgeschwindigkeit. (Um den 
richtigen Zusammenhang zwischen den Gréssen Qg, und Qy mit Q,, herzu- 
stellen, wurde aus den Resultaten von Abschnitt 53 ein provisorischer Wert 
fiir E herangezogen.) Nun ist man in der Lage, die Ausdriicke fiir Qx;, Oy 


und Q,, in Gleichung (5) einzusetzen und die Formel fiir den Wasseranteil des 
Eiswasseransatzes wie folgt anzugeben: 


W,, = 100 (1 OPT) xry B v1? + 0,785 Cy E (—tz) w, av  e6t 


4 (0,785 E wy d? vu — dm’/dx ) S 


(PT) xry sei dabei die Summe der beiden Druck-Temperatur-Funktionen. Wie 
spater gezeigt wird, ist dm’ vernachlassigbar; es ergibt sich damit der folgende 


. 
' 
) 
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vereinfachte Ausdruck 


Wr 100 (1 — AP IR ot) 


m — 


62,4 E wy; 2 (27) 


Es sei hier noch vermerkt, dass die Schritte, die zu Formel (26) fiihrten, fiir 
die vorgesehenen Experimente nicht unbedingt zu Recht bestehen; die An- 
nahmen tiber den Warmeiibergang und den Stoffaustausch beruhen darauf, 
dass auch die vereiste Kugel eine kugelige Umhiillende habe und deren Ober- 
flache glatt sei. Im Versuch diirfte dies kaum der Fall sein. Infolge der zulassigen 
Ahnlichkeitstheorie ist aber anzunehmen, dass die berechneten Verlaufe auch 


Q; 


L v 
5 10 15 20 m/s 


Figur 6 


Beitrage Qxz, Oy und Q,, zur Vernichtung der Gefrierwarme des angelagerten Wassers in Funktion 
der Luftgeschwindigkeit v; freier Wassergehalt wW,=5 g/m?, Kugeldurchmesser d = 2cm, Luft- 
temperatur 7, = — 10°C. 


den speziellen, von der Kugelform leicht abweichenden Verhaltnissen entspre- 
chen, wobei nur die Druck-Temperatur-Funktion (PT), mit einem zusatz- 
lichen Faktor @ korrigiert werden muss. Dessen Grésse kann jedoch nur expe- 
rimentell ermittelt werden. 

Gleichung (27) lasst eine Aussage iiber den maximal méglichen ieee 
anteil W,, zu, indem fiir w, > oo eine Grenzgerade gilt gemass 


Cw ee tr) 
lim Wy, ,99 = 100 (1 = ot) ' (28) 


Eine weitere Aussage ist zu erhalten fiir die wichtigen Bedingungen, wo 
W,, = 0 wird und der Eisansatz eine Temperatur von 0°C aufweist, das heisst 


HPT me Oe Ee l= 4) 


62,4 E wy i at a 


Die Lésungen von (29) nach ¢,, v, d und w, werden die Verhdltnisse charakte- 
risieren, bei denen eine Erhéhung der Temperatur, des Wassergehaltes, der Luft- 
geschwindigheit oder des Kugeldurchmessers die Ausbildung eines E1swasseran- 
satzes anstelle eines reinen Eisansatzes nach sich zient. 
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Auf die Darstellung der Formeln (27) und (29) sei an dieser Stelle verzichtet, 
sie wird im experimentellen Teil in den Figuren 11 bis 20 nachgeholt werden. 

Zusammenfassend kénnen die Hauptaussagen, zu denen die theoretischen 
Entwicklungen gefiihrt haben, wie folgt festgehalten werden: 

1. Der totale Warmeaustausch einer Eiskugel kann in seinen einzelnen 
Komponenten formelmassig erfasst werden. 

2. Bis auf einen konstanten Faktor #/F, der jedoch nicht stark von 1 ab- 
weichen kann ( und E liegen beide relativ nahe bei 1), ist es méglich, den 
Verlauf des Wasseranteils W,, in einem Eiswasseransatz in Funktion der Lutft- 
temperatur ¢,, des Kugeldurchmessers d, des freien Wassergehaltes w, der Luft 
und der Relativgeschwindigkeit Luft — Testkugel v wiederzugeben. 

3. Ferner ist die Asymptote bekannt, an die sich die Funktion W,, = [(w,) 
mehr oder weniger anschmiegen wird, je nachdem weniger oder mehr Tropfen 
von unserer Testkugel an die Umgebung abgegeben werden. 2 


3. Versuchsprogramm 


3.1 Allgemeines 


Auf experimentelle Weise sollen die Beziehungen zwischen den Ansatzge- — 


mengen an einer Testeiskugel und den sie produzierenden Vereisungsbedin- 
gungen ermittelt werden. Mit den gewonnenen Resultaten ist zudem die Giil- 
tigkeit der in Abschnitt 2 entwickelten, theoretischen Beziehungen nachzu- 
weisen. 

Um nicht von den Zufalligkeiten und den Launen der Natur abhangig zu 
sein, werden die notwendigen Versuche im speziell hierzu gebauten Hagel- 
versuchskanal ausgefiihrt. Diese Anlage erlaubt es, auch extremere Bedingun- 
gen zu erzeugen, als sie normalerweise in einer Gewitterwolke auftreten, und 
ergibt damit eine vorteilhafte Erweiterung des Variationsbereiches der Ver- 
suchsgréssen. 

Einerseits erzeugen und kontrollieren wir die Vereisungsbedingungen. Sie 
werden erfasst durch die Experimentparameter, andererseits untersuchen und 


charakterisieren wir das entstehende Produkt; die entsprechenden Aussagen 
werden durch die Messgrdssen gemacht. 


3.2 Die Experimentparameter 


Die Experimentparameter bezeichnen die Eigenschaften der das Messobjekt 
anstromenden Luft. Es sind dies vor allem die relative Geschwindigkeit Luft— 
Testobjekt, der Tropfenwassergehalt und die Tropfengréssenverteilung, die 
Lufttemperatur, der Luftdruck und der Eisbildungskerngehalt. 

Thre Variationsbereiche werden wie folgt angesetzt: Die relative Geschwin- 
digkeit Luft — Testobjekt oder, da das Testobjekt ortlich fixiert ist (siehe Ab- 
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schnitt 42), die Luftgeschwindigkeit v werde in drei Schritten von 6 m/s auf 
12, 18 und 24 m/s vergréssert. Damit werden auch kleinere und gréssere Ge- 
schwindigkeitswerte erfasst, als sie beim freien Fall auftreten. 

Der Tropfenwassergehalt w, sei verdndert zwischen 0 und 10 g/m? Luft 
oder mehr, wahrenddem die Tropfengréssenverteilung konstant bleibe. (Damit 
wird die Zahl der einzelnen Messungen wesentlich reduziert, ohne aber die Aus- 
sagemoglichkeit stark einzuschranken. Abgesehen von einer Veranderung der 
Einfangswahrscheinlichkeit diirfte beim Auftreten von Eiswassergemengen 
die Variation der Tropfengréssen keinen bedeutenden Effekt ergeben; anders 
ware es bei Vereisungsversuchen, wo alles Wasser gefrieren kann! Immerhin 
wird darauf geachtet, dass die mittlere Tropfengrésse von derjenigen in der 
Natur nicht allzuweit entfernt ist.) 

Die Lufttemperatur ¢; werde variiert in einem Bereiche von 0°C bis — 30°C. 
Diese Ausdehnung nach tieferen Temperaturen erweist sich jedoch nur sinnvoll, 
wenn der Gehalt der Luft an Eisbildungskernen, die in diesem Bereiche aktiv 
sein kénnen, praktisch gleich Null ist. 

Der Luftdruck # wird konstant gehalten und diirfte im Mittel 733 mbar 
sein, entsprechend der Lage des Forschungsinstitutes (2665 m ii. M.). 


3.3 Die Messgrossen 


Die Messgréssen sollen das Anfangs- und Endprodukt der Versuche, das 
heisst Ursprungspartikeln plus Eiswasseransatz, charakterisieren. Es sind dies 
der Durchmesser d und die Temperatur ¢, der Ursprungskugel sowie als 
Hauptmerkmal das Verhaltnis von Wasser zu Eis in einem resultierenden 
Ansatz. Dabei wird angenommen, dass die Gréssenénderung wahrend des 
Wachstums unserer Ursprungspartikeln und die damit verbundene Veranderung 
der Strémungsverhiltnisse dieses Messresultat in einem so geringen Masse ver- 
andern, dass die Abweichung noch im Bereiche der Fehlergrenze unserer 
Methode liegt. Dass diese Annahme zulassig ist, ergibt sich durch Versuche mit 
verschiedenen Kugelgréssen bei sonst gleichen Experimentparametern, wo nur 
eine relativ geringe Durchmesserabhangigkeit festgestellt werden kann. 


4. Messmethoden und Apparate 
4.1 Der Hagelversuchskanal 


Die Anlage, die das Erreichen und die Einhaltung der geforderten Werte 
der Experimentparameter gewdhrleistet, ist der Hagelversuchskanal. Dieser 


ist prinzipiell nichts anderes als ein im Messteil vertikaler, klimatisierbarer 


uw 


_ Windkanal mit geschlossenem Kreislauf. Seine Leistungsdaten wurden bereits 


frither publiziert [1] und werden daher nicht wiedererwahnt. 
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4.2 Der Einbau der Versuchsobjekte in die Mefstelle 


Da bis jetzt keine aerodynamischen Verhaltnisse geschaffen werden konn- 
ten, die ein einwandfreies Schweben von kugelahnlichen Partikeln gestatten, 
werden die Testobjekte der beschriebenen Versuche starr fixiert. Als Ursprungs- 
partikeln dienen kiinstlich gegossene Eiskugeln, in denen ein feiner Nylonfaden 
(nach Bedarf auch ein Thermoelement) eingefroren wird. Uber diesen Faden 


Figur 7 
Fixierung der Testeiskugel in der MeBstelle des Hagelversuchskanals (M 1: 1,33). 


ist ein feines Plexiglasréhrchen zu schieben, das an die Eiskugel angepresst 
werden kann. In dieser Position wird der Faden mit einem Stdpsel fixiert, 
wahrenddem das Réhrchen selber an einem Trager in der MeBstelle eingeklemmt 
wird (siehe Figur 7). Diese Art der Aufhangung bietet immerhin den Vorteil 
einer guten Isolation mit der Méglichkeit, auch Luftgeschwindigkeiten anzu- 
wenden, die von der freien Fallgeschwindigkeit des Testobjektes abweichen. 
Erste Versuche mit nur an einem Faden aufgehangten Objekten zeigten bereits 
bei 70% der Fallgeschwindigkeit gréssere, stérend wirkende horizontale Bewe- 
gungen und Pendelungen. Die Ausmerzung dieses Effektes ware wohl mit einer 
giinstigen Geschwindigkeitsverteilung iiber den Querschnitt der MeBstelle 
méglich, eine homogene Geschwindigkeit und ein genau definierter Feuchte- 
gehalt der an das Testobjekt angestrémten Luft sind aber vorzuziehen. 

An dieser Stelle muss auch ein Punkt erwahnt werden, der diese Versuche 
etwas vereinfacht: die Tatsache, dass der Eiswasseransatz aus thermodyna- 
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mischen Griinden eine Temperatur von 0°C aufweist. Das hat zur Folge, dass 
die Temperaturverhaltnisse geniigend- genau bekannt sind, sofern das Vorgehen 
beachtet wird, dass die Ursprungspartikeln vor der Einbringung in den Kanal 
eine Temperatur zwischen 0 und — 2° aufweisen. Eine merkliche zusatzliche, 
aber unerwiinschte Abfiihrung der Gefrierwarme kann damit nicht erfolgen. 


4.3 Die Messung der Windgeschwindigkeit 


Die Messung der Windgeschwindigkeit v an der MeBstelle geschieht tiber die 
Ermittlung der Druckdifferenz an den Enden des Effusors direkt vor der 
MeBstelle. Die Anwendung der Bernoullischen Gleichung ergibt dann die 
gewiinschte Aussage. Zusatzlich kGnnen auch Messungen mit einem Prandtl-Rohr 
ausgefiihrt werden — dieses Gerat ist jedoch infolge seiner Anfalligkeit gegen- 
iiber Vereisungen nur fiir Punktmessungen zu verwenden — im Gegensatz zur 
zuerst besprochenen Methode. 


4.4 Der Feuchtegehalt der Kanalluft, seine Erzeugung und Regulierung 


Die Erzeugung der Wassertropfen, die in den Kanal eingespritzt werden, 
geschieht mittels spezieller Farbspritzpistolen. Deren grosser Vorteil liegt 
darin, dass sie keine Pressluft bendtigen; die relativ feinen Tropfen werden 
bei hohen Drucken durch eine schwingende Membran produziert. Eine Gréssen- 
messung dieser Partikeln erfolgt durch Einfangen unterkihlten Nebels mit 
einem mit Paraffinél bestrichenen, — 30°C kalten Objekttrager (die auftref- 
fenden Tropfen gefrieren dabei im Moment des Auftreffens und kénnen nicht 
mehr koagulieren, sofern auch die mikroskopische Ausmessung bei negativen 
Temperaturen erfolgt) und zeigt, dass diese auch fiir Gewitterwolken als repra- 
sentativ angesehen werden diirfen (Durchmesser des mittleren Volumens 
@ ~504u). Je nach Wunsch kénnen 1, 2 oder auch 3 Tropfengeneratoren 
parallel geschaltet werden. 

Die Messung des freien Tropfenwassergehaltes geschieht tiber eine indirekte 
Methode, wobei die Tatsache ausschlaggebend ist, dass mit Ausnahme der nie- 
dersten Geschwindigkeitsstufe von 6 m/s iiberhaupt keine Vereisung der Kanal- 
wande stattfindet — im erwahnten Falle ist ein solcher Effekt zu erkennen; er 
ist jedoch bedeutungslos. - Wir versuchen daher, die total eingespritzte Was- 
sermenge durch Wagen zu ermitteln, wahrenddem weitere Messreihen die Ver- 
teilung des freien Wassergehaltes an der MeBstelle angeben sollen. Sie wird 


 gemessen, indem eine Sonde, unterteilt in 13 zylindrische Elemente von 2 cm 


Lange und 2cm Durchmesser nacheinander in den beiden Hauptrichtungen 
der quadratischen MeBstelle eingefiihrt wird. In der auszumessenden, vereisen- 


~ den Wolke wird nun die Gewichtsvermehrung der einzelnen Abschnitte durch 


_ Eisansatz bei — 20°C gewogen. Die resultierenden Verteilungskurven sind fiir 
- die verschiedenen beachteten Windgeschwindigkeiten in Figur 8 dargestellt. 
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Auf Grund der Kenntnis der absolut eingespritzten Wassermenge sowie deren 
Verteilung tiber den Querschnitt kann nun der absolute Betrag des freien 


Wassergehaltes angegeben werden. 


Wr 


Figur 8 


Mittlere Verteilung des freien Wassergehaltes iiber den mittleren Durchmesser dy der Mefistelle 
des Hagelversuchskanals bei verschiedenen Luftgeschwindigkeiten v. 


Die so ermittelten Eichkurven ergeben den Zusammenhang: total einge- 
spritztes Wasser — freier Wassergehalt am Ort des Testobjekts. - Wir miissen 
dabei jedoch eine unter Umstanden wesentliche Korrektur mit beriicksichtigen : 
der Feuchtegehalt der Kanalluft vor der Wassereinspritzung entspricht der 
Eissattigung bei der tiefsten Kanaltemperatur am Orte des Ammoniakver- 
dampfers oder Luftkiihlers (siehe auch [1]). An diesem Orte schlagt sich der 
letzte Rest der festen oder fliissigen H,O-Partikeln nieder, der tiberschiissige 
Wasserdampf kondensiert. Durch das eingespritzte Wasser wird nun die Luft 
vor der MeBstelle aufgewarmt und mit Wasserdampf gegeniiber den fliissigen 
Wassertropfen gesattigt. Die Distanz zwischen der Wassereinspritzung und der 
MeBstelle wurde nun so gewahlt, dass Zeit genug vorhanden ist, eine gegen- 
seitige Anpassung der Luft- und Tropfentemperaturen sowie des Wasserdampf- 
gehaltes zu gewahrleisten. Unter der Annahme resultierender stationarer Zu- 


stande ergibt sich eine Korrektur, die vom weiter oben angegebenen Resultat 
in Abzug zu bringen ist, von 


w=w—K, (30) 
K = Wey (ty) — Wsg(ty) - (31) 

w eingespritzte Wassermenge [kg/m3] 
K Korrekturwert zur Bestimmung des effektiven freien Wassergehaltes 


(kg/m?) 
Wsw Wasserdampfgehalt der Kanalluft beziiglich Wassersattigung [kg/m?] 
Wsr Wasserdampfgehalt der Kanalluft beziiglich Eissattigung [kg/m*] 
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tyy Temperatur an der MeBstelle [°C] 
ty Temperatur des Verdampfers [°C] 


Eine entsprechende graphische Darstellung ist Figur 9. Diese Zusammenhiange 
sind allgemein und gelten auch, wenn die Kanalluft zwischen Verdampfer und 
Wassereinspritzung aufgewarmt wird, sei es durch Heizung oder Warmelei- 
tung. 


K aes > = ae T q Se sf | 
gim*) ¢ 
| | | Ete 
2 aol 5 * 
rae a Sl 
-10 -20 -30 UG 
Figur 9 


Korrekturfaktor K zur Bestimmung des effektiven freien Wassergehaltes w, in Funktion der 
Verdampfertemperatur ty, bei verschiedenen MeBstellentemperaturen ty (= tz). 


4.5 Die Berechnung der Mefstellentemperatur 


Die Problematik der Temperaturmessung in einer unterkiihlten Wolke 
wird noch offensichtlicher, wenn diese nicht durch homogene Verhaltnisse ge- 
kennzeichnet ist. Da hilft auch keine integrierende Methode, wie es die Schall- 
geschwindigkeitsbestimmung eines reinen Tones und die damit erméglichte 


_ absolute Umrechnung auf die Temperatur darstellt. Die Vereisung jeder Art 
von Messfiihlern wiederum lasst ebenfalls keine geniigenden Werte zu. Aus 


: tscn 


Ate 


diesen Griinden wird die Errechnung der MeBstellentemperatur jedem andern 


Vorgehen vorgezogen. 
Die Gleichung fiir die Energiebilanz pro Kubikmeter Luft ergibt sich zu 


Wp Cy (ty — ty) = (ty — tscn) & Or + [Msw(tu) — @se(ty)] L- (32) 


In Worten heisst dies, dass die durch das pro m® Luft eingespritzte Wasser zu- 
gefiihrte Energie umgewandelt wird in eine Erwarmung der Luft, vermindert 
durch die Kalteproduktion beim Verdampfen von Wasser beim Angleichen an 
die neuen Sattigungsverhaltnisse. 

Die neu eingefiihrten Symbole bedeuten dabei: 
tw Temperatur des eingespritzten Wassers [°C] 
Temperatur, gemessen am Orte der ersten Umlenkschaufeln im 
Hagelversuchskanal (direkt vor der Wassereinspritzung) [°C] 


ZAMP .X1I/19 
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Cr spezifische Warme der Luft [kcal/kg] 
Or Luftdichte [kg/m*| 


Da die Kanalcharakteristik zeigt, dass die Schaufeltemperatur #5, gleich- 
gesetzt werden darf mit ¢,, der Verdampfertemperatur, so kénnen wir Glei- 
chung (32) trennen in eine Funktion der Verdampfertemperatur Q(t) und eine 
Funktion der MeBstellentemperatur Q(t), wobei Gleichung (32) tibergeht in 


Q(t) x Q (tyr) = 0, (33) 


Die Temperatur des eingespritzten Wassers sei dabei konstant. Diese beiden 
Funktionen sind dargestellt in Figur 10, wo auch die Differenz der gemessenen 


Alt) 
Qty) 
keal/h 


Bile 


tat, 
0 -10 -20 -30 ¢ 


Figur 10 
Darstellung der Temperaturfunktionen Q(t)) und Q(ty) fiir einen freien Wassergehalt w,= 10 g/m? 
und eine Temperatur des eingespritzten Wassers t,, = + 20°C. 


Verdampfertemperatur ¢, und der zu berechnenden MeBstellentemperatur ty, 
bei Q(t) = Q(ty) abgelesen werden kann. Figur 10 zeigt, dass bei héheren 
Temperaturen die zu beriicksichtigende Temperaturerhéhung gering ist — die 
durch das Wasser zugefiihrte Energie wird in erster Linie durch die Verdun- 
stung wettgemacht —, wahrenddem bei tieferen Temperaturen die Luft wesent- 
lich starker erwarmt wird. In diesen letzten Fallen ist der freie Wassergehalts- 
unterschied bei Erwarmung um einen Grad bedeutend kleiner, so dass der Ver- — 
dunstungsbeitrag abnimmt und die Warmezufuhr hauptsdchlich in eine Tem- 
peraturerhéhung umgewandelt wird. Bei einer Einspritzung von praktisch 0°C 
warmem Wasser wird die beobachtete Temperaturerhéhung wesentlich kleiner, 
als sie in Figur 10, der 20°C-warmes Wasser zugrunde liegt, dargestellt ist. 
Die berechneten Temperaturen scheinen damit gentigend gesichert zu sein. 
Aus dem Umstand, dass ein Temperaturwert mit einer absoluten Genauigkeit 
von + 1°C geniigt, ist auch eine Beriicksichtigung des Temperaturaustausches 


tiber den Querschnitt, hervorgerufen durch das Tropfenwassergehaltsprofil 
nicht vonnoten. 
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4.6 Die Bestimmung des Wasseranteils von Eiswasseransatzen 


Die Bestimmung des Wasseranteils in Eiswasseransaitzen kann am besten 
tiber eine kalorimetrische Messung der Enthalpie gegeniiber einer Temperatur 
hdher als 0°C vorgenommen werden. Zu diesem Zwecke tauchen wir das 
vereiste Testobjekt, das Endprodukt unserer Versuche, in ein Dewar-Gefiss, das 
teilweise mit einer Sperrfliissigkeit angefiillt ist. Die resultierende Temperatur- 
anderung dieses Bades erlaubt es, auf den urspriinglichen Wassergehalt der 
vereisten Partikel zu schliessen. Unter der Annahme, dass diese eine Temperatur 
von 0°C aufgewiesen habe, ergibt sich die Aussage, dass die durch das vereiste 
Testobjekt zugefiithrte Energie E; gleich sein muss wie die Reduktion der 
innern Energie der Sperrfliissigkeit E.,. Diese beiden Gréssen lassen sich aus- 
driicken durch 


Ep = Gg S + Gp Oy tay; (34) 

und 

Esp = (tsp.o — tui) Gsp Esp » (35) 
wobei 
Gy Gesamtgewicht des im Testobjekt enthaltenen Eises [kg] 
Gr Totalgewicht des Testobjektes [kg] 
tui Mischtemperatur [°C] 
tsp-o Ursprungstemperatur der Sperrfliissigkeit [°C] 
Gsp Gewicht der Sperrfliissigkeit [kg] 
Csp spezifische Warme der Sperrfliissigkeit [kcal/kg°] 


Die Kombination der Gleichungen (34) und (35) ergibt 


G Gy Sp a 
— E= (ts, . ty) > > a Kp. (36) 


Kp Korrekturfaktor, charakteristisch fiir das verwendete Dewar-Gefass 
(wird durch Eichung ermittelt) 


Da wir jedoch am Wasseranteil W,, des Ansatzes interessiert sind und nicht 


-am Wasseranteil der gesamten Testpartikeln, muss die folgende Operation 


beriicksichtigt werden: 


y a Gw _ Gr GE 37 
Wy AO0 reat = LO) so SG (37) 

mv, Wasseranteil des Eiswasseransatzes [%] 

eG A Gewicht des Eiswasseransatzes an die urspriingliche Eiskugel [kg] 

Gx Gewicht der urspriinglichen Eiskugel [kg] 

2G, Gewicht des fliissigen Wassers [kg] 


4, 54 


Da bei dieser Methode die Warmekapazitat und die Tragheit in der zeit- 
lichen Anpassung des Dewar-Gefasses eine merkbare Rolle spielen kénnten, wur- 
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den Eichkurven aufgenommen, die den Wert und die zeitliche Abhangigkeit 
des Korrekturfaktors Kp ergeben, mit dem das scheinbare Gewicht des Eises 
G;, multipliziert werden muss, um den korrekten Betrag G,; zu erhalten. Die 
Differenz aus Laboratoriumstemperatur und Ausgangstemperatur der Sperr- 
fliissigkeit wird in allen Versuchen praktisch gleich gewahlt und die jeweiligen 
Phasen der Vorbereitung und der Messung gleich gestaltet. Es sei auch noch 
erwahnt, dass als Sperrfliissigkeit stets Wasser verwendet wird; zur Messung 
der Temperaturen und Temperaturdifferenzen dienen eine Kompensations- 
messbriicke und ein Spiegelgalvanometer. Alle Gewichtsbestimmungen werden 
auf zwei verschiedenen Mettler-Prazisionsanalysenwaagen ausgefiihrt. 

Die Voraussetzungen zu diesem oben beschriebenen Ermittlungsverfahren 
scheinen auf den ersten Blick an eine homogene Temperatur von 0°C der 
Ursprungskugel auch nach der Vereisung gebunden zu sein. Effektiv findet 
aber auf deren Leeseite eine Abkiihlung statt. Durch Einbau von Thermo- 
elementen an der Riickseite und im Zentrum verschiedener Testkugeln konnte 
jedoch nachgewiesen werden, dass dieser Effekt ganz unbedeutend ist. (Die 
értlichen Warmeiibergange kénnen auch nach den von FROssLING [10] ermit- 
telten Warmeiibergangszahlen abgeschatzt werden. Dabei ist jedoch zu beriick- 
sichtigen, dass sich die Grenzschicht durch Aufnahme von Gefrierwarme er-— 
warmt und dadurch der Warmeiibergang auf der Leeseite des Testobjektes 
reduziert wird!) Dies bedeutet auch, dass frei fallende, rotierende Eiskugeln in 
erster Nadherung gleiche Wasseranteile W,, in einem Eiswasseransatz auf- 
weisen wie gleich grosse, jedoch aufgehangte Testkugeln, sofern die Experi- 
mentparameter dieselben sind! 

Zur Bestimmung des Wasseranteils von Eiswassergemengen wird noch eine 
zweite Methode angewendet: Das Endprodukt eines einzelnen Vereisungsexperi- 
mentes wird in eine Zentrifuge gebracht und ausgeschleudert. Dabei muss 
darauf geachtet werden, dass dies bei 0°C geschieht und keine weiteren Phasen- 
libergange das Resultat verfalschen. Der Wasseranteil eines Eiswassergemenges 
kann auf diese Art rein durch Wagung bestimmt werden, was gegeniiber der 
kalorimetrischen Methode eine wesentliche Vereinfachung der Experimente 
bedeutet. Allerdings bleibt eine Unsicherheit bestehen, da unter Umstdnden in 
geschlossenen Kavernen vorkommendes, fliissiges Wasser in einer Zentrifuge 
nicht entfernt werden kann. Dies ist mithin der Grund, der dazu fithrte, beide 
Methoden parallel nebeneinander zu verwenden. 


5. Messungen und Resultate 
5.1 Vorgehen 


Die Auslegung des ganzen Messprogramms ist abhangig von den Experi- 
mentparametern. Wahrend es die Anlage, der Hagelversuchskanal, erlaubt, die 
Temperatur und die Windgeschwindigkeit auf gewiinschte Werte einzustellen, 
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ist es nicht méglich, den freien Wassergehalt, der durch die Spritzpistolen 
eingefithrt wird, beliebig zu wahlen. Kleinste Finfliisse auf die Diisen konnten 
die Durchflussmenge verandern, auch darf zur Erlangung einer homogenen und 
stets gleichbleibenden Tropfengréssenverteilung keine Verstellung an den 
Pistolen vorgenommen werden. (Es hat sich gezeigt, dass eine Verianderung 
der Spritzmenge stets verbunden ist mit einer Verinderung der Durchmesser 
der erzeugten Tropfen.) Zusatzlich verandert sich bei konstanter Einspritzung 
der Wassergehalt indirekt proportional mit der Luftgeschwindigkeit. 

Das sind die Griinde, die eine Verwendung von drei Spritzpistolen nétig 
machten, die, einzeln oder wahlweise kombiniert, verschiedene freie Wasser- 
gehalte erzeugen kénnen. Bei niederen Luftgeschwindigkeiten ist jedoch ein 
Einsetzen aller Tropfengeneratoren nicht mehr méglich, da damit ein zu hoher 
freier Wassergehalt verbunden ware; die daraus resultierenden Eiswasser- 
gemenge kénnten nicht mehr geniigend genau bestimmt werden. 

Infolge der geschilderten Umstande wird bei einer bestimmten Lufttem- 
peratur und einer bestimmten Kanalgeschwindigkeit sowie einem ebenfalls 
konstanten Testkugeldurchmesser der Wasseranteil des Eiswasseransatzes in 
Funktion des freien Wassergehaltes bestimmt. Nach Ausfithrung einer solchen 
Messreihe wird die Luftgeschwindigkeit auf einen neuen Wert eingestellt und 
eine neue Reihe begonnen. Sind alle vier gewiinschten Geschwindigkeiten 
durchgemessen, so kann fiir eine neue MeBserie eine andere Temperatur gewahlt 
werden. Die sich aus diesem Vorgehen ergebenden Resultate sind in Abschnitt 
(5.2) zusammengestellt. Auf eine Variation des Durchmessers der urspriing- 

lichen Testeiskugel wird generell verzichtet, also nur Einzelmessungen iiber 
die Abhangigkeit des Wasseranteils von demselben ausgefiihrt (Abschnitt 5.3). 


Das Vorgehen, das bei einer Einzelmessung (Bestimmung des Wasseranteils 
_ kalorimetrisch!) eingehalten wird, geschieht chronologisch wie folgt: 
, a) Giessen der Testkugel mit dem gewiinschten Durchmesser, 
b) Einstellung des Hagelversuchskanals auf eine bestimmte Temperatur und 
; eine gegebene Luftgeschwindigkeit, 
c) Vorbereitung der Sperrfliissigkeit zur kalorimetrischen Bestimmung des 
Wassergehaltes der vereisten Probe und Temperaturbestimmung, 
d) Wagen der Sperrfliissigkeit, 
~ e) Vorbereitung der Testkugel (Bestimmung des Durchmessers, des Gewichtes 
sowie Kontrolle der Temperatur, Aufbringung der Distanzhiille und Fixie- 
rung des Nylonfadens, siehe Abschnitt 4.2), 
f) Ablesen von Kanaldaten (Verdampfertemperatur, Temperatur vor der 
Wassereinspritzung, Temperatur an der Mefstelle, Luftgeschwindigkeit), 
___ eventuelles Nachjustieren einzelner Werte, 
_ g) Einbau der Testkugel in die MeBstelle des Hagelversuchskanals, 
__h) gleichzeitige Inbetriebsetzung der Wassereinspritzung, 
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i) unmittelbar nach dem Abschalten der Spritzpistolen sofortiger Ausbau des 
vereisten Objektes, 
k) Entfernung der Distanzhiilse, 
1) Wagen des vereisten Testobjektes, 
m) Einbringen desselben in die Sperrfliissigkeit des Kalorimeters, 
n) Abstoppen der Schmelzzeit (zu Kontrollzwecken), 
) Durchmischung der Sperrfliissigkeit plus Testobjektwasser inklusive Be- 
stimmung der Mischtemperatur, 
p) Bestimmung der in den Kanal eingebrachten totalen Feuchte, 
q) Nachkontrolle der Kanalbedingungen. 


O 


Das Vorgehen bei der Einzelmessung andert sich, wenn der Wasseranteil 
mit der Zentrifuge und durch Wagung bestimmt wird. Die Punkte c und d 
fallen weg, wahrenddem die Positionen m,n und o ersetzt werden durch das 
Auszentrifugieren der vereisten Testpartikeln und die nachfolgende Wagung. 
Die kritischste Stelle des Versuchsablaufes ist die Zeitdauer, welche der Ausbau 
der vereisten Kugel aus dem Kanal, ihre Wagung und ihre Uberfiihrung in das 
Kalorimeter oder die Zentrifuge bendtigen. Mit einiger Routine kann dieser 
Vorgang sehr beschleunigt werden, so dass normalerweise keine Wasserverluste 
durch Austropfen aus dem Eiswassergemisch auftreten. Die Warmezufuhr in 
der Zeit, da das Objekt mit der Laboratoriumsluft in Kontakt ist, diirfte ver- 
nachlassigbar sein, héchstens vermindert sie die Unterkiithlung der Ursprungs- 
partikeln (ca. 2°C). 


5.2 Messresultate 


Die auf Grund des in Abschnitt 51 beschriebenen Vorgehens ermittelten 
Ergebnisse sind in den Tabellen 1 und 2 eingetragen und in den Figuren 11, 12 
und 13 dargestellt. Tabelle 1 enthalt dabei die Messwerte, wie sie mit der kalori- 
metrischen Methode bestimmt wurden, wahrenddem Tabelle 2 die Wasseran- 
teile, bestimmt durch Zentrifugierung und Wagung, enthalt. Es zeigt sich dabei: 
a) Wird der prozentuale Wasseranteil W,, des Eiswasseransatzes in Funktion 

des freien Wassergehaltes aufgetragen, so lasst sich bei konstanten Werten 

der Geschwindigkeit und der Temperatur eine Verbindungskurve durch die 

Messpunkte legen; diese Kurve bestimmt zugleich einen freien Wasserge- 

halt, bei dem der Wasseranteil 0% betragt ; 

b) fiir wachsendes w, wird ein Grenzwert lim W,, angestrebt, der kleiner ist 
als 100% ; 

c) limW,, geht mit abnehmender Temperatur nach kleineren Werten iiber; 

d) der Wassergehalt W,, nimmt allgemein fiir kleinere Geschwindigkeiten 
kleinere Werte an; 


e) der Wasseranteil des Eiswasseransatzes sinkt allgemein bei abnehmender 
Temperatur ¢,. 


He oe « 


avi 


~ AM 
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Tabelle 1 


Messwerte des Wasseranteils W,,, und deren Abweichungen AW vom theoretischen Wert 

(in Prozenten des totalen Ansatzes G4) in-Funktion des freien Wassergehaltes w,, bei ver- 

schiedenen Luftgeschwindigkeiten v-und Temperaturen ¢;, Kugeldurchmesser d — 2 cm. 
Bestimmung des Wasseranteils gemdss kalorimetrischey Methode. 


tr —10°C —20°C —30°C 
» wy | Wm | AW | wy | Way | pees ag w,, | AW 
mis | gim’ | % |%Ga | sim? | % | %Ga | gm? | o% | %G4 
24 | 18] 10 1+ 2 4° 16 0 o| 23 ae | 
4,5 | 46 Onl 6%] 038 (no 5l.90.0 |) 24 Jo 4 
8,9 | 65 Ont 2.5 0 Gilt 16,35) B4de elsbe <3 
SBA Sar 5, ok wh 7420 la) Beas el 4s alg 
3.4 | 32 |-12 417.6 | 64 ‘Ila 2 
a Oe Uhl tee) rhe hee gat 0.0. bad ibaa, 8 
oo theas |= acb-6.6 igo 2b gz lessee 4a 
ae Pesord}at3 «(| 25,6. 25 lor 6.54 weoti-se ele y 
Be a6H eG. 4 44.95 ire oheens Brogan ha ae Sete 
B41 20 Del ig a 4 seule sg 
BAL bese be Ff PRS O is 
i4,8 | 25 |—10 
OO ENT SG aes ae ee a) ed ey TN 
DT db sGh i 2) G4 levee i 8.1) 143 on hen 3 
SWUM 1S: fs 4-1) 15.61% 34 or 6:71.35 ated 
on) ay 11290 \) 25 [22 12's 5 
0} 46. 14 
6 | 9,7 | Gk eee | 12.0 34. | o | 15,9 | 2 Vol 
3.6 heen 9 oes ae 8 Loos 33 ees 
el Ng le a I Oe 9 Og ae ee gs 
10,0 | 51 |- 9 


Die absoluten Werte lassen sich fiir alle Temperaturen relativ gut durch 
Kurven darstellen, wobei zu bedenken ist, dass einerseits die Messung des 
freien Wassergehaltes w, um 15% falsch sein kann, anderseits der Wasseranteil 
W,, gemass Abschnitt 46 bisweilen mit einem absoluten Fehler von bis zu 15% 
des totalen Ansatzes G, behaftet ist. Speziell interessant ist dabei die Uberein- 
stimmung der Resultate, die einerseits durch kalorimetrische Messung, ander- 
seits durch Zentrifugierung und Wagung gewonnen wurden! Dadurch werden 
die Ergebnisse wesentlich gefestigt, wenn auch beiden Methoden die Eigenschaft 
zugeschrieben werden kann, eher etwas zu kleine Werte zu ergeben. 

Dass gerade bei der héchsten Temperatur von — 10°C die gréssten Streu- 
ungen der Messwerte auftreten, ist nicht verwunderlich. Das angelagerte 
Wasser gefriert hier sehr langsam und hat dadurch Zeit, an den Ort des Krafte- 
gleichgewichtes (bestimmt durch den aerodynamisch bedingten statischen 


- Luftdruck, den hydrostatischen Druck, die Schwerkraft und die Oberflachen- 
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Tabelle 2 ' 


Messwerte des Wasseranteils W,, und deren Abweichungen AW vom theoretischen Wert 

(in Prozenten des totalen Ansatzes G4) in Funktion des freien Wassergehaltes w,, bei ver- 

schiedenen Luftgeschwindigkeiten v und Temperaturen 77, Kugeldurchmesser d = 2 cm. 
Bestimmung des Wassevanteils mit Zentrifuge und Wagung. 


tr —10°C —20°C are 
v Wy Wi, | AW Ws | Won | AW wy Wm | AW 
m/s | g/m®)| % | %Ga | e/ms | % | %Ga | g/m? | ee MGA 
24 BE. | eeSGenl aad 1,0 0 0 6,1 4 ae 
3,5 | 45 0 64 12s TSG 4,4 1 tee 
B35 | GO.n es 8,3 | 43 0 een ip rake 0 
ii | BO el =e G O83 50 8 eed 
2504p ees 55d |) ASD Nes 
18 | 6.9.11, 067 ila 2.5) For een eones 7,0-| > 10 sleeiass 
5-0 t P 52° Tle | a Gol aes eel 6,6 7 a al 
13,5] 72 “|or 45 | 14,5¢4) ge TN een es 
9.6 |: 63 |X 7,5).11,0 1S i 2,5 
12° | 10;7°) "62" |=: 8,5) F,0 aieeaee 2 eto ate 0 
Bo bs Si 0 H13,3a eae eee 7.8 SOME | ARLO 
793933 | 
6°) 12:5 | 60 ‘le 4.) tomas ees | 13.4 16 ai 7 
8,5| 45 110" 4 06 Se eeeornl ee 
S45) 4B Wetng SiR al eee a ae 
10:0 |. 59 -|="%2 


Wr 
% 
24m/s 
— 18m/s 
> 1am/s 
6m/s 
50 


e 2% m/s 
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Figur 11 


Experimentelle Werte und theoretisch-empirische Kurven iiber den Zusammenhang Wasseranteil _ 


W,,, des Ansatzes — freier Wassergehalt w, der Luft, bei verschiedenen Luftgeschwindigkeiten v 
und der konstanten Lufttemperatuy ty = —10°C, Kugeldurchmesser d = 2 cm. 


ie a » « 
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Figur 12 
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Experimentelle Werte und theoretisch-empirische Kurven iiber den Zusammenhang Wasseranteil 
W,, des Ansatzes — freier Wassergehalt w, der Luft, bei verschiedenen Luftgeschwindigkeiten v 


und der konstanten Lufttemperatur tr = — 20°C, Kugeldurchmesser d = 2 cm. 
a © Amis 
% 
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Figur 13 


10 


Experimentelle Werte und theoretisch-empirische Kurven tiber den Zusammenhang Wasseranteil 
W,,, des Ansatzes — freier Wassergehalt w, der Luft, bei verschiedenen Luftgeschwindigkeiten v 


und der konstanten Lujttemperatur tr = — 30°C, Kugeldurchmesser d = 2 cm. 


spannung) zu fliessen. Bei hohen Luftgeschwindigkeiten gelangt es an den Ort 
des tiefsten statischen Druckes an der Kugelperipherie und bildet dort einen 
eigentlichen Ring; wahrend dem bei kleineren Geschwindigkeitswerten die 
Schwerkraft fiir die Wasseransammlung im Staupunkt auf der Luvseite der 
_Luftstrémung massgebend ist. (Die bereits friiher erwahnte Erscheinung [11], 
dass fiir die Form des resultierenden Ansatzes in erster Linie die Vereisungs- 
bedingungen verantwortlich sind, tritt hier besonders deutlich zutage!) Dies 
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hat zur Folge, dass in Funktion der Geschwindigkeit eine starke Variation der 
sich bildenden Verformung — ausgehend von der urspriinglichen Kugel — auf- 
tritt. Ihrem wechselnden Einfluss auf die Einfangswahrscheinlichkeit kann nur 
begegnet werden, indem relativ kleine Ansatze erzeugt werden, wodurch jedoch 
der zugehérige Messfehler erhoht wird. 

Punkte, die auf die Abszissenachse zu liegen kommen (siehe Figuren 11, 12 
und 13), sagen aus, dass alles Wasser bereits wahrend des Wachstums gefroren 
ist. Sie konnen nicht weiter ausgewertet werden, da es sich um reine Eisansatze 
handelt, fiir welche die theoretischen Uberlegungen gemiss Abschnitt 2 nicht 
mehr geniigen. 

Es sei hier auch darauf hingewiesen, dass es nicht ohne weiteres selbstver- 
standlich ist, wenn auch die Messungen bei ¢; = — 30°C in unser Bild passen, 
sind doch die Resultate eindeutig an die Bedingung gekniipft, dass die an die 
Testkugel angelagerten Wassertropfen die Umgebungstemperatur angenommen 
haben, ohne zu gefrieren. Rechnerisch konnte gezeigt werden, dass eine Abktih- 
lung der an die Testkugel angeschwebten Wassertrépfchen auf die Lufttempe- 
ratur in der zur Verfiigung stehenden Zeit und bei den in Betracht zu ziehenden 
Durchmessern gewahrleistet ist. Dies ist auch experimentell dadurch nachge- 
wiesen worden, dass anstelle des im Normalfall eingespritzten 0°C-warmen 
Wassers 20°C-warmes Wasser eingespritzt wurde und kein merkbarer Einfluss 
auf den Wasseranteil der so produzierten Eiswasseransatze nachgewiesen 
werden konnte. Zusatzlich ergab sich dabei, dass die eingespritzten Wasser- 
tropfen nicht nur die Lufttemperatur praktisch annehmen, sondern auch ihren 
unterkiihlten Zustand beibehalten und nicht gefrieren. Dies ist einerseits den 
im Hagelversuchskanal eingebauten elektrostatischen Filtern [12] zu verdan- 
ken — ein Ausschalten der Filter lasst die Wassertropfen 6fters gefrieren, indem 
vom Verdampfer mitgerissene Eispartikeln gefrierprozessauslésend wirken -, 
anderseits dem Umstand zuzuschreiben, dass das in den Kanal eingespritzte, 
normalerweise unbehandelte Wasser keine Fremdsubstanzen enthalt, die bei 
den auftretenden Luftgeschwindigkeiten in der kurzen Strecke zwischen den 
Wasserdiisen und der MeBstelle eisbildend wirken. 


5.3 Anpassung der Theorie an die Messresultate 


Die Anpassung der in Abschnitt 2 entwickelten Theorie wird wesentlich er- 
leichtert durch die Tatsache, dass ein Ablésen fliissiger Wassertropfen aus det 
Wasserhaut eines vereisten, kugeligen Testobjektes nicht beobachtet werden 
konnte. Diese Aussage ist aber fiir Wasseranteile W,,, die grésser sind, als sie 
im Experiment vorkommen, nicht gesichert, auf alle Fille spielt aber diese 
Erscheinung niemals die Rolle, die ihy in der Literatur oft zugeschrieben wird. 


Selbstverstandlich miissen Beobachtungen an schmelzenden Eispartikeln 
davon ausgenommen werden. 


 -— a 
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Dies bedeutet, dass wir ohne Bedenken Gleichung (27) verwenden diirfen, 
wobei nur noch die Aufgabe verbleibt, den Koeffizienten 0/E auf Grund von 
Messungen zu bestimmen. Hierzu sei der Punkt W,, = 43%, w,= 8,3 g/m, 
v = 24 m/s, t; = —20°C (Figur 12) gewahlt, der eingesetzt im Ausdruck (27) 
fiir 0/E den Wert von 1,48 ergibt. Die damit errechneten Kurven sind in den 
Figuren 11, 12 und 13, in denen die Messpunkte einzeln aufgefiihrt wurden, als 
ausgezogene Linien eingezeichnet; ebenfalls eingetragen sind auch die berech- 
neten Grenzwerte fiir lim W,,. Der Vergleich zwischen Messpunkten und be- 
rechneten Kurven zeigt, dass die Ubereinstimmung unter Beriicksichtigung der 
verschiedenen Fehlerméglichkeiten, gemass den in Tabellen 1 und 2 angefiihr- 
ten absoluten Abweichungen, als sehr befriedigend bezeichnet werden kann. 
0/E hangt prinzipiell von allen Parametern (Partikeldurchmesser d, Geschwin- 
digkeit v, Temperatur ¢,, freier Wassergehalt w,) ab; um so auffalliger ist es, 
dass ein einziger konstanter Wert zur Charakterisierung der Messresultate 
geniigt. (Die Tendenz, dass die Messwerte mit sinkender Temperatur die 
theoretisch ermittelten Kurven erreichen oder sogar tiberschreiten, ist nicht 
deutlich und gesichert genug, um daraus eine Abhangigkeit von #/E abzu- 
leiten.) 

Die Formel, die die Ergebnisse beschreibt, erhalt mit der Einfiihrung des 
gefundenen Korrekturfaktors #/E fiir einen Luftdruck von 733 mbar (ent- 
sprechend der Lage des Institutes) die folgende Form: 


W,, = 100 — 2,37 (PT) py 4-2 v8 w! + 1,26 t, . (38) 


Dieses Resultat soll uns jedoch nicht davon abhalten, die effektiv méglichen 
Einfliisse auf #/E kurz zu diskutieren. Bestimmend fiir 3/E sind speziell der 
Raumwinkel der Ansatzflache auf der Ursprungskugel inklusive die Form des 
Ansatzes (die ja stets von der fiir die Rechnung zugrunde gelegten, idealen 
_ Kugelbegrenzung in Form und Rauhigkeit abweicht), die Einfangswahrschein- 
lichkeit und die jeweilige Abhangigkeit dieser Gréssen von den Vereisungsbe- 
dingungen. Eine Basierung auf die Ergebnisse von Messungen anderer Autoren, 
wie zum Beispiel LANGMUIR [13], ist nicht unbedingt ratsam, da alle bis jetzt 
gemachten Untersuchungen iiber Einfangswahrscheinlchkeiten darauf beru- 
hen, dass das angelagerte Wasser direkt am Ort der Anlagerung vereist wird. 
Bei den Eiswasseransitzen ist jedoch der Umstand, dass das angelagerte 
Wasser sich auf der Oberflache in Funktion der Umgebungsbedingungen ver- 
schieben kann, entscheidend. Solche komplexe Wachstumstendenzen sind bis 
jetzt nicht in Erwagung gezogen worden, auch ist die dazugehdrige Unter- 
suchungsmethodik noch nicht geniigend entwickelt. 
Eine experimentelle Ermittlung der Durchmesserabhangigkeit des Wasser- 
anteils eines Einsatzes zeigt einen eher betonteren von d abhangigen Verlauf 
; (W,, ~ a), als er durch Formel (38) gegeben ist (W,, ~ d~/*), unbertihrt des 
- Umstandes, dass die Messungen innerhalb der angegebenen Fehlergrenzen 
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[+ 15%] mit der Theorie iibereinstimmen. Diese Tendenz ist verstandlich mit 
Riicksicht auf die Abhangigkeit des Faktors #/E vom effektiven Durchmesser d. 

Die Ermittlung der Druckabhangigkeit der Formeln (27) oder (38) scheitert 
an der Komplexitat des Faktors #/E, wogegen alle andern Gréssen ohne 
weiteres auf beliebige Drucke umrechenbar waren. Um diese Schwierigkeiten 
zu meistern, miisste vorerst eine Kennzahl bestimmt werden, die bei Vereisun- 
gen gleiche Tropfentrajektorien miteinander in Beziehung bringt und damit 
analog ist zu der Grésse, die RusH und WARDLAw [14] benutzen. Nur sollte die 
neue Kennzahl fiir Eiswassergemenge gelten, bei denen die an ein Ursprungs- 
objekt angelagerten Tropfen nicht am Orte ihrer Anlagerung gefrieren. Trotz- 
dem darf aber in erster Annaherung angenommen werden, dass bei einer 
Absenkung des Druckes bis auf 200 mbar der Koeffizient #/E nicht allzu stark 
variieren wird. Da die Druck-Temperatur-Funktion (PT),;y nur schwach vom 
Druck abhangt (Figur 3), so diirfte fiir eine erste Abschatzung Formel (38) 
geniigen. Bei der Umrechnung auf frei fallende Partikeln wird sich aber zusatz- 
lich ein anderes Verhaltnis zwischen v und d ergeben, da bei konstantem 
Durchmesser die Fallgeschwindigkeit bei kleiner werdenden Drucken nach 
grossern Werten strebt. 


Auf Grund der Messungen und den durch die Theorie erméglichten Extra- _ 


polationen sollen abschliessend all die Abhangigkeiten des Wasseranteiles in 
einem Eiswasseransatz zusammengestellt werden. Infolge der relativen Un- 
iibersichtlichkeit der Gleichung (38) sei W,, in Funktion aller vier Experiment- 
parameter dargestellt. Figur 14 zeigt den Wasseranteil in Funktion des freien 
Wassergehaltes w, bei einer Luftgeschwindigkeit von 24 m/s und bei verschie- 
denen Temperaturen, wahrenddem Figur 15 die Zusammenhange darstellt in 
Funktion der Lufttemperatur und bei verschiedenen Wassergehalten. Figur 16 
zeigt W,, als Funktion der Geschwindigkeit v fiir verschiedene Temperaturen 
und Wassergehalte. Wahrend in den eben erwahnten Figuren 14, 15 und 16 
der Durchmesser der Testkugel konstant gehalten wurde, ist er in Figur 17 


r 


Wal 
% 


60} 


5 10 15 g/m® 
Figur 14 
Darstellung des Wasseranteils W,, in Funktion des freien Wassergehaltes Ws, Parameter Lufttemi- 
peratur 77; Luftgeschwindigkeit v = 24 m/s, Kugeldurchmesser d = 2 cm. 
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Figur 15 
Wasseranteil W,, in Funktion der Lufttemperatur ¢z, Parameter freier Wassergehalt w,; Luftge- 


' schwindigkeit v = 24 m/s, Kugeldurchmesser d = 2 cm; W, = Grenzwert des Wasseranteils. 


Figur 16 
Wasseranteil W,, in Funktion der Luftgeschwindigkeit v, bei verschiedenen Werten'der Lufttempe- 
ratur tp und des freien Wassergehaltes Wy, Kugeldurchmesser d = 2.cm. 


Figur 17 


- Wasseranteil W,, in Funktion des Durchmessers der vereisten Kugel, berechnet fiir freien Fall der- 
_selben, bei verschiedenen freien Wassergehalten w, und Lufttemperaturen tz; Luftdruck 
P p = 733 mbar. 
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variiert fiir verschiedene Werte der Lufttemperatur und des freien Wasser- 
gehaltes. Alle Angaben fiir Durchmesser unterhalb 1,0 cm und oberhalb 2,5 cm 
stellen dabei Extrapolationen dar. 

Fiir die Erérterung von Figur 17 wurde es als vorteilhaft erachtet, den 
Durchmesser mit jener Grenzfallgeschwindigkeit in Zusammenhang zu bringen, 
die sich bei frei fallenden Eiskugeln einstellen wiirde. Die dafiir zu machenden, 
zusatzlichen Annahmen sind die folgenden: Hagelkorndichte oy = 900 kg/m? 
(in Ubereinstimmung mit gemachten Messungen [15]), Luftwiderstandskoeffi- 
zient c,, = 0,5. Auf diese Weise errechnete Wasseranteile gelten aber nur unter 
dem Vorbehalt allfalliger Korrekturen von seiten der Strémungsdynamik [11]. 


Generell steht damit fest, dass abnehmende Werte von Temperatur, frecem Was- 
sergehalt w,, Lufigeschwindigkeit v und Kugeldurchmesser d stets eine V erminderung 
des Wasseranteils W,, zur Folge haben. 


Um die praktische Bedeutung der Eiswassergemenge bei atmospharischen 
Vereisungen noch besser hervorzuheben, seien die Vereisungsbedingungen zu- 
sammengestellt, bei denen gerade noch Eisansatze entstehen, die keine Wasser- 
einschliisse mehr aufweisen. Zu diesem Zwecke sei die Gleichung (29) verwendet, 


woraus sich die in Figur 18 dargestellten Verhaltnisse ergeben. Es zeigt sich dabei, — 


Fis -Wasser 


-20 -30 
Figur 18 
iets er Sepa” fiir verschiedene Luftgeschwindigkeiten v, die das Gebiet der reinen 
iswassergemengen in Funktion des freien Wassergehaltes w, und der Luft- 
temperatur ty trennen, Eiskugeldurchmesser d = 2 cm. 4 


°C 
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dass die Grenzkurven fiir verschiedene Luftgeschwindigkeiten die Koordinaten- 
ebene (w,,¢,) in zwei getrennte Gebiete aufteilen. Im obern Gebiet herrschen 
Bedingungen, die zu Eiswassergemengen fiihren, im untern Gebiet jene, die 
reine Eisansatze bewirken, wahrenddem die Grenzkurven selber den geometri- 
schen Ort darstellen fiir reine Eisansitze, deren Temperatur 0°C betragt. 
Der Eiskugeldurchmesser wurde hierbei mit 2 cm als konstant angenommen. 

Figur 19 zeigt-dieselben Grenzkurven, jedoch fiir den Fall, wo der Kugel- 
durchmesser variiert, die Geschwindigkeit aber konstant bleibt. 

Schliesslich charakterisiert Figur 20 jene Grenzfalle, die sich fiir frei fallende 
Partikeln ergeben, unter den gleichen Bedingungen, wie sie der Figur 17 zu- 
grunde liegen. Figur 20 kann auch herangezogen werden, um vergleichende 
Betrachtungen beziiglich der Hagelbildung anzustellen, indem sie Grenzkurven 
fiir verschiedene Kugeldurchmesser frei fallender Eispartikeln zeigt. Kleine 
Eiskérner bendétigen bet jeder unter dem Gefrierpunkt liegenden Temperatur weit 
grossere freie Wassergehalte als grosse Hagelkorner, damit Eiswassergemenge 
anwachsen kénnen. Diese Erscheinung riihrt speziell vom spezifisch gréssern 
Stoff- und Warmeaustausch her, welcher den kleinern Kugeln eigen ist, zudem 
wird er durch die Abhangigkeit der Grenzfallgeschwindigkeit von der Quadrat- 
wurzel des Partikeldurchmessers verstarkt. 


Figur 19 


Darstellung der Grenzkurven fiir verschiedene Eiskugeldurchmesser d, die das Gebiet der reinen 
_ Bisansatze von den Eiswassergemengen in Funktion des freien Wassergehaltes wy und der Luft- 
f temperatur tz trennen; Luftgeschwindigkeit v = 18 m/s. 
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6. Schlussbemerkungen 


Die Messungen iiber das Wachstum von Eiswassergemengen an Eiskugeln 
bei atmosphdrischen Vereisungsbedingungen zeigen die grosse Bedeutung 
dieser bis jetzt unbekannten Wachstumsméglichkeit von Hydrometeoren, die | 


g/m 


Fis-Wasser 


i he | eke 
-10 -20 -30 ep 
Figur 20 


Darstellung der Grenzkurven fiir frei fallende, kugelige Partikeln verschiedenen Durchmessers d, 
in Abhangigkeit des freien Wassergehaltes w,und der Lufttemperatur tz ; Luftdruck = 733 mbar. 


auch bei den haufigsten schweren Flugzeugvereisungen eine massgebende Rolle 
spielen kann. Speziell zeigt Figur 20, dass primar in jeder Wolke, in der Ver- 
eisungen durch Anlagerung von unterkiihlten Wassertropfen stattfinden, Eis- 
wassergemenge bildende Bedingungen notwendig vorhanden sein miissen. Erst 
die genaue Kenntnis dieser Vorgange wird mithelfen, eine brauchbare, den 
atmospharischen Bedingungen dquivalente Hageltheorie zu schaffen und die 
physikalische Seite der Hagelbildung befriedigend abzuklaren. 

Es ist meine angenehme Pflicht, Herrn Prof. Dr. R. SANGER und Herrn 
Prof. Dr. G. Buscu fiir ihre férdernde Anteilnahme an dieser Arbeit meinen 
herzlichsten Dank auszusprechen; Dank gebiihrt auch dem Direktor des 
Schneeforschungsinstitutes, Herrn Dr. M. DE QUERVAIN sowie Herrn Dr. Tu. 
ZiNGG fiir ihr stetes Interesse und die Unterstiitzung der beschriebenen Unter- 


, 
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Suchungen. Meinen beiden Mitarbeitern, Herrn P. RAss und Herrn J. VON 
NIEDERHAUSERN bin ich fiir ihre Mitwirkung bei der Ausfiihrung der Messungen 
ebenfalls verpflichtet. Die Arbeit wurde im Rahmen des Forschungsprogramms 
der Eidgenéssischen Kommission zum Studium der Hagelbildung und der 
Hagelabwehr an deren Forschungsstelle Weissfluhjoch am Eidgenéssischen 
Institut fiir Schnee- und Lawinenforschung ausgefiihrt ; es standen hierzu auch 
wesentliche Kredite des Schweizerischen Nationalfonds zur Verfiigung. 
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Summary 


The first experiments with the Swiss Hail Tunnel [1] have shown that the 
growth of hailstones under atmospheric icing conditions does not occur in the way 
usually assumed. It is quite normal to observe still liquid water in the resulting 
ice deposit. Figure 1 shows a model of this new observed mechanism while Figure 2 
shows an ice framework without the liquid water normally contained in it. 

These first findings led to the experiments and the theory which are described 
and developed in this paper. The final goal was to determine the water content of 
the growing ice-water mixtures as a function of the generating icing conditions and 
the characteristics of the primary ice particles. For the calculations and the experi- 
ments we used in both cases ice spheres with an initial temperature of 0°C. The 
theoretical total exchange of heat is given in formula (26) or — as the experiments 
have shown that no water normally escapes from the iced spheres — by formula (27). 
The results of the experiments are given in Figures 11, 12 and 13, which also indicate 
the calculated values. To achieve a correlation between theory and experiment, it was 
necessary to determine an unknown value #/E where E means the collection effi- 
ciency and # a corrective factor. #/E was found to be constant at all different 
values of the free water content, the air velocity, the air temperature and the 
diameter of the test particles (ranging from 1:0 to 2:5cm). The accuracy of the 
measurements is -- 15% of the total deposit; the mean volume diametre of the 
water droplets, accrued on the test particles, is 50. 
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Summarizing it can be said, that increasing values for temperature, free water content 
air velocity, and the diameter of the ice spheres lead to increasing values for the liquid 


water component in the resulting deposit. ; ; 
Of particular interest are the conditions which determine the borderlines be- 


tween deposits of dry ice and deposits of ice-water mixtures. These are indicated in 
Figure 20 which gives the corresponding growing conditions as functions of the 
free water content and the air temperature. Figure 20 is valid for freely falling ice 
spheres of different diameters and shows, that the deposits consisting of 1ce-water 
mixtures ave growing preferably at higher temperatures and on greater particles. In 
-addition we can also recognize the large range of atmospheric conditions which 
produce aqueous hailstones. These results should demonstrate the great importance 
of the thermodynamics of such ice-water mixtures as icing deposits on freely 
falling particles. 


(Eingegangen: 6. Dezember 1959.) 


Steady-State Heat Conduction in a Circular Cone’) 


By Roxuro Mukti and ELI STERNBERG, Providence?), R.I., U.S.A. 


1. Statement of Problem. Solution in Complex Integral Form 


Consider a circular cone of semi-opening angle 9) (0 < 6, < a). Let rectan- 
gular Cartesian coordinates (x, y, z) be chosen in such a way that the origin 
is located at the vertex and the z-axis is coincident with the axis of the cone. Let 
(r, 8, y) be spherical coordinates defined by the mapping 


x=rsinOcosp, y=rsinOsing, z=rcos6, 


1) 
(Sri co, VSO Sa, Oe ans 
and introduce the auxiliary notation 
p=cos0, py=cos,. (2) 


In the rotationally symmetric problem to be treated presently the tempera- 
ture I’ is independent of the longitude gy and obeys the boundary conditions 


Tr, p) =T, for OSrSa, Tir,p)=0 for a<r<oco, (3) 
while throughout the interior of the body the temperature ‘field T(r, p) must 


1) The results communicated in this paper were obtained in the course of an investigation 
conducted under Contracts Nonr 562(20) and Nonr 562(25) of Brown University with the Office 
of Naval Research in Washington, D.C. 


*) Division of Applied Mathematics, Brown University. 
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satisfy LAPLACE’s equation, which governs the steady-state heat conduction of 
a thermally homogeneous and en solid, Thus 
(sid & 2 
“Or2- oe x Y ae te a , 5s male a p ) a | — (4) 


must hold in the region 0 <7 < oo, 0<4@< @. Finally, we adjoin the 
regularity requirements 


T = o(r-), Le =o(r*) as roo, | 
ar i (9) 
Fie th), a O(r) as r>0. | 


The Dirichlet problem characterized by (3), (4), (5) is most conveniently 
attacked by means of the Mellin transform?). Suppose 


ce [TC p) 1 dr, (6) 


so that T(p, s) designates the Mellin transform with respect to 7 of T(r, ), 
s being the transform parameter. Integration by parts, in view of (5), now yields 
the identities 


[i nar=—s is [Beran ser tes, (7) 
provided the real part of s obeys the restriction 0 < Re(s) < 1. On multiplying 
(4) by 7? and applying the transform (6) one obtains, by virtue of (7), the 
Legendre equation 
d 
ap 


[a9 37) + se- Feo. 8) 


~ On the other hand, the boundary conditions (3) in the transform domain take 


ae 


on the form 


F by, s\ et (9) 


The solution of (8) that meets (9) and corresponds to a temperature field 
which is bounded in the physical domain under consideration, is evidently 


given by 
: Tya® P,_(p) 
Le ee ao) 


where P_, is the Legendre function of the first kind and of degree s — 1. 


3) See Dortscu [1] and SnEppoNn [2] for general discussions of this transform. Numbers in 


brackets refer to the list of publications at the end ot the paper (page 314), 
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According to the inversion theorem for the Mellin transform we have 
t T —s 
TV, 2) = 357 / T(p, s) r-* ds , (11) 


and substitution from (10) into (11) yields 


ra a | i ACA (12) 


The formal solution (12) for the desired temperature distribution depends 
upon the choice of the line of integration Re(s) = c in the complex s-plane. 
Clearly, the solution is not affected by the particular choice of c as long as the 
position of the line of integration is varied within one and the same strip of 
regularity of the integrand in (12). The selection of the appropriate strip, in 
turn, is dictated by conditions (5). 

With a view toward establishing the admissible values of c, we observe 
first that the integrand in (12) is a meromorphic function of s, the poles of which 
coincide with the zeros of the function 


F(s, po) = $ Pa(Po) - (13) 


It is known‘) that all zeros of P.(p)) are real and simple for every choice of 
py (— 1 < py < 1). Also, from the Mehler-Dirichlet integral representation of 
Legendre functions’), 


0 
PGi = 2p 
0 


cosa — cos 


(0<6<2), (14) 


whence P,_,(f9) is positive for 0<s <1. Consequently, the region 
0 < Re(s) < 1 is contained wholly in a strip of regularity of the integrand in 
(12), the actual width of this strip depending on the specific value of f,. If we 
take c in the range 0 <<c <1, s conforms to the restrictive assumption under- 
lying (7). It remains to be shown that this determination of the parameter c 


is also sufficient for (5) to hold; this is indeed the case as will become apparent - 


in the succeeding section. 


2. Residue Expansion. Solution in Real Integral Form 


We now evaluate the definite integral (12) with the aid of the theorem on 
residues and thus deduce an infinite series solution of the problem at hand. 


To this end we assume c fixed in the interval 0 < c <1 and complete the 


4) Macnus ard OBERHETTINGER (3], p. 91. 
5) See WHITTAKER and Watson [4], p. 815. 
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contour of integration in the half-planes which lie to the left or the right of 
Re(s) = c, according as OS r<a or a<r< ow. Specifically, the two closed 
contours selected, say C and C’, are defined as follows: they have in common the 
straight-line segment Re(s) = c, — @ S< Im (s) < @ (e > 0); in addition, C and 
C’ consist of the respective semi-circles |s — c| = g, Re(s) <c and |s —c| =, 
Re(s) 2 c. Evidently, the integration around C must proceed in a counter- 
clockwise sense, while C’ is to be traversed in the opposite direction. 
It follows from the asymptotic expansion 


(0) Vso Ft foelle S)0— 4] + te 


| (15)°) 


as |s]| oo («S0Sn-¢€, ¢>0) 


in which J" is the gamma function, that P._,(p)/P,_1(f)) remains bounded as 
| s | oo. Bearing in mind (12), as well as the definitions of C and C’, we 
therefore conclude that 


T(r, #) =lim = f f(s) ds (0<r<a), | 
: (16) 
T(r, b) =lim 55; fils) ds (a<r<ov), | 
@ 
where 
fo = BEM. om 


The poles of the meromorphic function f(s) are coincident with the zeros of 
F(s, po) defined by (13), all of which, we recall, are real and simple, and lie 
outside the interval 0 < s <1. These zeros, with the exception of the one at 
the origin, in turn, coincide with the zeros of P_,(po). Let the positive and 
negative zeros of P(p,) be designated by s,(po) and sy(po) (k = 1,2, .--), 


respectively, and be ordered in accordance with 


gee 6 se a SS (18) 

_ From the recursion formula 
Pb) = P..(6) (19) 

we then have the relations 
ge sp 1 (Se ee (20) 


ERDELY! [5], p. 162. Recall also that P,(1) = 1. 
Here r, p, and py) may be regarded as parameters. 


; ; 
a 
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On applying the residue theorem to the contour integrals in (16) one 
obtains, in view of (18), (19), (20), the expansions 


rena ni+s ae ey! (| esr<an 
} (21) 


Ss, +1 


Tp) = Todt nampa (Zeon, dacine cee | 


In order to demonstrate the absolute and uniform convergence of the two 
infinite series appearing in (21) for the corresponding range of y and for 
0 <0 <4, it is evidently sufficient to show that 


eons fee! 


= 
olay Cea ee (22) 


Si, 
provided 0 < 6 < 6) and 0 < 6) < a. To see that (22) holds one merely needs 
to recall that P(p) is bounded for 0 < 0 < z if s is real, bear in mind that 


dP(p,)/ds cannot vanish at a simple zero s = s,(po) of P.(p 9), and make use 
of the asymptotic expansions 


aP(po) _ V2 Oo sin[(s + 1/2) + 3 2/4] +0(s3) (s>0) 
ds Vx s sin 0, (23)8) 


as $s —> co. (eS 0; — 6, -e— 0); 


and 


sy(bo) = G- (k — 4) — +00) 
leis (24) 
as k>oo («S65 a-—€, €> 0). 


Since the functions 


P.(p) rt, P.(6) 7 (25) 
are both harmonic, the same is true of the individual terms in each of the 
series (21). Further, it is readily verified that the series obtained from (21) 
by a termwise application of the Laplacian operator (referred to spherical 
coordinates) remain uniformly convergent in the respective domains 0 < 7 < a, 
0S6< 6, anda<r< 0,0 <6< 6). Consequently T(r, p) in (21) satisfies 
LAPLACE’S equation (4). According to the definition of Sup) (k= 172, 37am 
8) This estimate follows from (14) by differentiation with respect to s and subsequent recourse : 


to an asymptotic result given by ERDELYI [6], p. 49. 
PAL (71: 


Yo. ae es 


ee 
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P, (bo) = 0, whence (21) conform also to the boundary conditions (3). Finally, 


(18) and the convergence properties of (21) permit the conclusion that the 
temperature field under examination meets the regularity conditions (5), as 
anticipated at the end of Section 1. This completes the verification of the in- 
finite series solution (21) and, by implication, justifies the formal complex 
integral representation (12) if0<c <1. 


A real integral representation for the desired temperature distribution may 
be obtained from (12) by taking c = 1/2 and introducing the change of variables 


S ae = +7¢&. (26) 
In this manner there results after some elementary manipulation 
re yfa of. 
rer.) = 22% [eter 0) a, (27) 
0 


where 


Y : Y 
cos (¢ log) — 2 €sin (e log) Po vy2sie(?) 
g@ 1/4 P_tjosig(Po) 


sé, 7, p) = (28) 


The Legendre functions occurring in (28) are real-valued, as is apparent from 
(14). It is not difficult to show without recourse to the previous convergence 
considerations that this alternative representation of the solution fulfills the 
requirements (3), (4), (5). Such an independent validation of (27) is, however, 
superfluous. Equation (27), in contrast to the individual infinite series in (21), 
holds for the full range 0 < 7 < ow. 


3. Evaluation of Solution. Numerical Results 


In this section we discuss procedures for the numerical evaluation of the 
results derived in the preceding section. Turning first to the series solution (21), 
we observe that its evaluation necessitates a knowledge of the zeros s,( 9) of 
Pp); in addition, one requires the values of P,(p) and dP,(pp)/ds at s = s,(p). 

A scheme for the computation of the roots of P.(f9) = 9 (regarded as an 
equation in s) was devised by MacponaLp [9]. Pat [7], [8] developed a more 
efficient method of computation and actually determined s,(9) for 
k = 2, 3, 4, 5, 6 and several choices of fy. The family of roots corresponding to 
k = 1, overlooked by Pat, was supplied later on by Horton [10]. For the sake 
of convenience, we list in Table 1 relevant values of s,(p9) drawn from [7], [8], 
[10]. Additional entries in this table may be computed readily by Pat's 
procedure. 
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Table 1 
Values of s;,(cos 9p) 

9% | m/12 1/6 m/4 n/2 

iP 8:45 4-13 Asie) 1-00 

2 20-60 10-03 6.52 3-00 

3 32-55 16-02 10-51 5-00 

4. 44-53 22:01 14-51 7:00 

| 5 56°53 28-01 18-50 9-00 
6 68:52 34-01 22:50 11-00 


The calculation of P, (p) on the basis of the Mehler-Dirichlet integral (14) 


is awkward because the integrand is singular at 1 = 0. It is more expedient to 
employ instead the following proper integral representation which may be 
deduced from (14) through integration by parts. 


0 
P(p) = nd he [cos (cosA — cos 6)1? 
7 Sin 
6 “) (29) 0) 
~ - (s + 2) (s — 1) (cosA — cos 6) 2] cos |(s ++ 3) al ah. 
An analogous formula for dP.(f9)/ds is reached directly from (29) by differentia- 
tion with respect to s and need not be recorded here. 

For the evaluation of the solution in the integral form (27), (28) the values 
of the Legendre function P_,).,;<(p) are needed. Suitable tables of this function 
do not appear to be available. Equation (29), however, yields the computa- 
tionally convenient integral representation 


sin? @ 
6 


6 
P_1)2.:e(P) ee 42° [cos (cosA — cos 6)1?2 | 
(30) 

ae S (= a 3) (cosA — cos 6)s| cosh (E A) da, | 


which holds for0 << @9<2z,0<&< ow. In particular, for € = 0, (30) implies 
the well-known result?) 


P_spld) = = K(sin 3), (31) 


where K is the complete elliptic integral of the first kind. — 


10) This equation is valid for 0 < 9 < 2. Recall that Pea = 1: 
11) Cf. Erpéiyi1 [5], p. 174, 


ey Ew 
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The numerical evaluation of the proper definite integral (30) for fixed 
6 (0 < 6 < a) presents no difficulties-Once the requisite values of P_ 12+ie(D) 
are known one may proceed to the evaluation of the improper integral (27) 
by numerical means. In this connection it is helpful to have an asymptotic 
estimate of the corresponding truncation error. Such an estimate may be 
reached with the aid of the asymptotic formula 


Pe dt ef + OF] as E00 (e SP Sa—e, > 0), BY 
: V2 &sin0 


which follows from (15). Indeed (28), (32), and integration by parts, lead to 


[et ae =a te | 
[ 


sin 6 
& 


[ (4 — 4) sin (Slog =) + log cos(s, log ~)] exp[— fo (09 — 8)] 


x 


Eo | (0. — 8)* + (log) ] 


+ Of{€5% exp[—&) (0) — 8)]} as £00 (O<1r<oo, eS0<6,e>0),] 


ee ; Y Y tA 
2/2 sin 6, [% sin (é, log =) ae log— cos (é log" )| exp (— 


V5 [5% + (log) | 


[e714 = sae 
+ 0 [é>3? exp(—&,0))] as E200 (0<1r<oo). J 
Figure 1 shows the variation of the temperature along the axis (0 = 0) of 
the cone for 0, = 2/6, i.e. for an opening angle of sixty degrees. These numerical 
values were obtained on the basis of the solution in the real integral form (27), 
(28) by the scheme just outlined. The underlying numerical integrations were 
carried out on an I. B. M. 650 electronic computer. The same figure also 
displays the dependence of T upon y at 6 = 0 for 6) = 7/2, i.e. for the special 
case of a half-space which is exposed to a uniform temperature Tj over a 
circular subregion of its plane boundary, the remainder of the boundary being 
held at zero temperature. A closed solution in terms of elliptic integrals to the 
latter problem was given previously in [11]. The results found in [11] reduce 
to an elementary form on the axis of symmetry. Introducing the special values 
pp = 0, P =1 into the residue expansion (21), and making use of known prop- 
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0,=2/6 


= 
0 0:5 1-0 1-5 2:0 25 
r/a —_—_> 
Figure 1 
Variation of temperature along axis of cone. 


erties of the Legendre functions of the first kind, one arrives at 


¥ 


Tr WLT l! = Ta (0<1<co) (35) 


when 6, = 2/2, which is in agreement with [11]. 


Acknowledgment 


The authors are greatly indebted to Mrs. L. PAuL who conducted the 
numerical work involved in this investigation. 


REFERENCES 


[1] G. Dorrscu, Handbuch der Laplace-Transformation, Vol. 1, Birkhauser, Basel 
1950, 
[2] 1. N. SNEDDoNn, Fourier Transforms, McGraw-Hill, New York 1951. 
[3] W. Macnus and F. OBERHETTINGER, Formeln und Sdtze fiw die speziellen 
Funktionen dey mathematischen Physik, Springer, Berlin 1948. 
[4] a T. WuItTTaKeER and G. N. Watson, A Course of Modern Analysis, Cambridge 
946. 


[5] Higher Transcendental Functions, edited by A. Erp&xyt, Vol. 1, McGraw-Hill, 
New York 1953. 


[6] A. ERDELYI, Asymptotic Expansions, Dover, New York, 1956. 


Vol. XI, 1960 Kurze Mitteilungen — Brief Reports — Communications bréves Se 


[7] B. Pat, On the numerical calculation of the roots of the equations en) 0 anc 
d/du Pj"(u) = 0 regarded as equations in n. Bull. Calcutta Math. Soc., 9, 2, 85 
(1917). 

[8] B. Par, On the numerical calculation of the roots of the equations P™(u) = 0 and 
d/du Pi"(u) = 0 regarded as equations inn [Part II}. Bull. Calcutta Math. Soc. 70, 
3, 187 (1918). 

[9] H. M. Macponatp, Zeros of the spherical harmonic P,™(u) considered as a function 
of nm. Proc. London Math. Soc., 37, 264 (1899). 

[10] C, W. Horton, Note on the zeros of P,™ (cos @) and dP™ (cos 0)/d0 considered as 
functions of n. Bull. American Math. Soc., 53, 2, 153 (1947). 

[11] E. STERNBERG and E. L. McDowELt, On the Steady-state Thermoelastic Problem 
for the Half-space. Quart. Appl. Math., 74, 4, 381 (1957). 


Zusammenfassung 


Diese Arbeit befasst sich mit der Bestimmung des stationaren Temperaturfeldes, 
das in einem Kreiskegel durch eine unstetige Verteilung der vorgegebenen Ober- 
flachentemperatur hervorgerufen wird. Die Temperatur der Mantelflache wird 
oberhalb und unterhalb einer festen Entfernung von der Kegelspitze konstant an- 
genommen. Dieses Problem wird mit Hilfe der Mellin-Transformation streng 
gelést. Die so erhaltene Lésung wird in einer spater erscheinenden Arbeit auf die 
Bestimmung der zugehérigen Temperaturspannungen angewendet. 
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Note on the Effect of a Centre of Dilatation in a Semi-Infinite 
Elastic Solid with its Plane Face Fixed 


By Smritt Kana BHOWMICK, Calcutta, India*) 


Introduction 


The object of this note is to find the effect of a centre of dilatation in the form ofa 
nucleus of elastic strain (say) inside a semi-infinite elastic solid having its plane 
face rigidly fixed. The method of solution is simple, and the results are obtained in 
a closed form. 


Solution 


We take the plane surface of the isotropic elastic solid as the plane z = 0 with 
- the origin O on the surface. The axis of z is drawn into the body. If there be a 
centre of dilatation at A (0, 0, c), we have the displacement (uw, v,, w,) due to it in 


1) Jadavpur University, Department of Mathematics. 
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the form (LovE, p. 127) 


0 i) 0 1 
[Uy, V1, wy | = P | Ties ) (1) 


Ox’ Oy? Oz 
where P is a constant depending on the strength of the nucleus and elastic constants, 
while 
Ri = #74 ye (2 yer (2) 

Taking a nucleus of dilatation of equal and opposite strength at the image point 

B (0, 0, — c), we have the corresponding displacement (ug, vg, W,) given by 
0 0 0 1 

[Us, V9, W| a P ae: oy ? rea vee, ? 


(3) 


where 


R= w+ 9? + (2 +0). (4) 


Combining the two systems of displacement (1) and (3) we find that on the 


surface z = 0 
LEPC 
Un Uy NO, Utes =n Oecd @ +O = pag» (5) 
in which ie 
R? = (4? + y? + 6?) . 


Since by hypothesis, the plane boundary is rigidly fixed, we are to find the ~ 


third system of displacement (uw, v, w) which will satisfy the equation of equilibrium 
and nullify the displacements (5) on the surface z = 0. 


8 (0 0, -C) 


Centre of dilatation in an elastic solid. 
The equations of equilibrium in the absence of body forces are 


V [u, Uv, w]) a ai an 2 o [5 4 Oy x | ¢ , 4 
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where A is the dilatation and o is Potsson’s ratio. Putting 


: OF 
A p=2 (1 — 20) —— (7) 
such that 7? F = 0, we obtain from (6) 
0 0 0 
72 [y — (72 Fa — 
P* tt, 0, 2) =f see oP Mae Ee 
Then we can put 
[u, v, w} = [25 Cine gf F 8 
7%) = Ox?” Oy? 4 al + [Pi Pa Ps) » (8) 


where (,, 2, ~3 are harmonic functions having prescribed values on the surface z = 0. 
Since we are to nullify the results given in (5), we can take in this case gy, = 0, 
~, = 0. Then from the relation 


Ou ov . Ow , OF 
Bg ee Gyo gg Laas 
we shall have 
al 
latinas 3 9 
tp (Sao) ie 


Assuming the harmonic function 


0 Sg 2P (z+) 
ge us ) 10 
Pe 2P > (R,] Fg eae ~*~ 


we find that the displacements (5) on the surface z= 0 are cancelled. Hence from 
(7) we get 


| ES eErG 


, 11 
3—4o Vie Ae) 


Substituting these values of y, and F in (8) we get (uv, v, w). The resultant displace- 


_ ments (U, V, W) are therefore given by 


AAPA 


1 1 62 (2+ ¢) } 
U = uy + uy tu= Px | R + Re rte tae 
il 1 62 (2+ am 
Van+H+0=Py| Re T Re (3 = 40) R3 , (12) 
z2—C # ace 22 1 
W=u,+u,+w=P ere a ere 
o2 | | 
(3= 40) R§ 
a : 
~ On calculating stress components from U, V, W, it is found that on the surface 
faa 
Ge = Vo 0), (13) 
while } ; 
4EP(l-—o ee | if 
22],-0 = (7 +6) (3 — 40) a= R 0 
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Zusammenfassung 


Es wird das Verschiebungsfeld ermittelt, das im elastischen Halbraum mit 
fixierter Grenzflache durch ein Dilatationszentrum erzeugt wird. 


(Received: February 8, 1960.) 


Permanent Replicas of the Crystalline Structure of Hailstones 


By Pace R. ParnTerR}), and VINCENT J. SCHAEFER’), Schenectady, N.Y., U.S.A. 


The internal structure of hailstones has attracted the attention of natural 
scientists for many years. The recognition of the layering which often occurs dates 
back at least to Vota [1]*), who, in 1800, made a sketch of the onionskin-like 
appearance of stones which fell at Como, Italy. The more recent work in this field 
is tipified by the elegant cross section of hailstones collected in Switzerland published — 
in color by List [2]. Due to the evanescent nature of hailstones, very little of a 
quantitative nature has thus far been attempted in the United States, although 
there are few areas in the world where such storms are more frequent or the storms 
more violent and damaging. 


Developing Interest in Hailstorms in America 


With the advent of cloud modification and evidence that profound changes can 
be initiated in cumulus clouds when properly seeded with dry ice and silver iodide 
{3, 4], a number of serious studies of hailstorms have been conducted in Nebraska, 
Colorado, New York, West Virginia, Oregon, California in the United States, in 
Alberta, Canada, in the vicinity of Lake Garda in Italy, the vineyard country of 
France, and in Switzerland where an excellent, comprehensive research study has 
been underway for more than five years under the auspices of the Federal Commis- 
sion for the Study of Hail Formation. The initial work has consisted almost entirely 
in the operation of silver iodide ground generators. The rather casual efforts of 
earlier years has more recently been directed toward seeding with intensive net- 
ey of ground generators, groups of aircraft and in some cases self destroying | 
rockets. 

In much of the work of the past, nothing but very qualitative studies have been 
made of the end product of the hailstorm — the hailstone. This is understandable 
since the random nature of severe storms rarely permits the interested scientist to 
be present during or immediately after the storm when observations are best made. 


This has been particularly true in the intense storms of the Great Plains of America 
where the stones often fall in isolated areas. 


1) Davis High School, Davis, California. 
2\ 


) Rensselaer Polytechnic Institute, Troy, New York. 
8) Numbers in brackets refer to References, page 326. 
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A new approach to this problem is now possible and its feasibility was demon- 
strated conclusively during the Summer of 1959, 


Research Projects on Hailstorms at the Loomis School 


In the Spring of 1959, the senior author was asked to plan and direct a summer 
school for talented high school students sponsored jointly by the American Meteoro- 
logical Society and the Loomis School of Windsor, Connecticut, with support from 
the National Science Foundation. 

In drawing up a series of research projects to be pursued by the students, several 
research activities were directed toward studies of the hailstorms which occasionally 
ruin the tobacco crop of the mid-Connecticut area. These included the observation 
of storms by radar, the collection and preservation of hailstones, and the develop- 
ment of methods for preserving their internal structure by a replica technique. All 
three projects were worked on by students during the summer. 


Studies of the Replication of the Structure of Hailstones 


In order to develop the replica technique for preserving the crystalline structure 
of hailstones, hydromites were formed in a cold chamber. These were made by 
placing water in a polyethylene bag suspended above a cold chamber operating 
at —20°C and prepared so as to have a very slow leak, permitting one drop to fall 
about every half minute. The drops froze and gradually built up a cylinder of ice. 
By varying the rate of falling droplets, it is possible to control the freezing rate at 
the surface of the hydromite thus producing various sizes and shapes of crystals 
within the ice. These formations were used in developing techniques for cutting, 
smoothing, polishing, and etching ice surfaces. After suitable etch figures developed, 
the surface was coated with a cold (—5°C) 3 percent solution of polyvinyl formal 
dissolved in ethylene dichloride after the manner described by SCHAEFER [5, 6]. 


Acquisition of Hailstone Samples 


In order to assure that hailstones were available for the critical evaluation of 
the techniques developed, aid in obtaining hailstones was sought from T. HENDER- 
son of Weather Modification Company who was working on a hailstorm prevention 
project in northeastern Colorado. In mid-July an air express shipment of fairly 
large (2—2!/, cm diameter) hailstones arrived at Windsor. Subsequently, a second 
shipment of considerably larger stones (up to 5 cm diameter) was obtained from 
Grand Island, Nebraska, through the excellent cooperation of Mrs. GAYLORD 
WALLACE of Burwell, Nebraska and S. FLAHERTY of the Grand Island office of the 
U.S. Weather Bureau. 


Technique for Preparing Hailstone Replicas 


After making a number of experiments, the authors recommend the following 
technique for the preparation of hailstone cross section replicas: 

1. A block of plastic foam of the styrofoam type is prepared having the dimen- 
sions 3cm xX 3cm X 5cm. 

2. A depression is made in one surface of the foam about 2 cm in diameter and 
0-5 cm deep. Three symmetrically placed holes are poked in the foam at the bottom 
~ edge of the depression to anchor the ice which is to be used for adhesive. 
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3. About one half milliliter (ml) of water is placed in the depression in the 
styrofoam. The amount of water used depends on the size of the hailstone. 

4. The cold hailstone (about —20°C) is placed in contact with the water, causing 
it to freeze immediately. So long as the assembly is kept colder than 0°C, the stone 
can thenceforth be handled and manipulated without danger of breakage. The 
plastic foam is such a good insulator that if the stone is normally kept at about 
—20°C there is no danger of structural change. 

5. The hailstone may be cut in half with a jeweller’s saw, an ordinary saw, a 
mitre saw, or a band saw, depending on the nature and size of the stone. The 
localized melting which occurs seems to have no effect on the structural features 
of the stone, provided sawing is done while the ice is at a temperature of —20°C. 

6. After sectioning, the surface of the stone is smoothed, using 2/0—100 and 
6/0-220 grade sandpaper‘). These are prepared as strips 6” x 22” long. They are 
clamped on a flat board with the working surfaces about 20” long. As they become 


loaded with ground ice new pieces are substituted. Care must be taken that deep — 


scratches from foreign objects do not occur during the grinding procedure, other- 
wise the replica will not be of good quality. After smoothing the surface is polished. 

7. Polishing may be done with chamois, paper felt, or velvetine. This is mounted 
adjacent to the sandpaper for convenience. While it is desirable to accomplish the 
polishing without causing surface melting, it has been observed that when an ice 
surface is polished so vigorously as to cause localized melting, the etch figures 
develop without apparent modification. This presumably is due to the presence of 


the crystalline structure directly below which immediately replicates the underlying ~ 


crystal patterns even if a thin melt layer forms momentarily. 
8, After polishing the surface, the hailstone is held at a temperature of about 
—20°C for 15 to 30 minutes in a cold chamber or a cold room at a humidity lower 


than ice saturation. Localized evaporation from the ice crystal matrix forming the ~ 


surface occurs as ice molecules sublime away, leaving behind crystalline patterns 
identical in visual appearance to the figures which appear on etched metallographic 
specimens. 

9. After the etch figures become visible on the surface of the hailstone, a layer 
of replica solution should be flowed over the surface. This should consist of 5% 
polyvinyl formal dissolved in ethylene dichloride. This solution should have a tem- 
perature of —5°C, and, after coating, the surface should be inclined at 45° to 
drain off excess material. The specimen should be placed in a chamber with some 
ventilation and held at a temperature of about —20°C until the solvent has 
completely evaporated. This normally occurs in about 20 minutes. 

10. The replica film may be recovered in several ways. By holding the hailstone 
at —5°C for a day or so, sufficient sublimation occurs through the film to loosen 
the replica film from the surface. This film may be carefully lifted and transferred 
to a glass slide for mounting. In the event additional sections are to be made of the 
stone, a cut can be made about 2 millimeters below the film on the band saw and 
the thin section melted from the replica film by placing it on a glass slide with the 
film-coated surface against the glass. If the stone is to be sacrificed, the same 
procedure may be followed, or it may be placed in a tray of water with the film on 
the surface, and, after the stone is melted, the replica slid onto a glass slide for 
mounting. Another effective procedure is to keep the stone in a warm room until 
surface melting occurs, at which time the film may be slid off the surface onto a 
glass slide wet with water. The replica film is preferably mounted with the smooth 
upper surface in contact with the glass plate. 

*) 


We use Behr-Manning Adalox No. 6931 and No. 364. 


— 
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Figure 1 


Hailstones from two severe storms in the Great Plains of the United States: 


Upper: Grand Island, Nebraska, June 27, 1959. 
Lower: Sterling, Colorado, July 17, 1959, 


ZAMP X1/21 
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Figure 3. Detail shown by replica at center of hailstone. Replica of etched surface varies in 
brightness as replica is rotated. 


a THAW vale HNC} TU INPUT HNN b 


MM 11 2 3 4 5 


Figure 4a. Appearance of replica in normal light. Figure 4b. Appearance of replica between 
crossed polaroids with air inclusions showing polarized refraction. 
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11. After drying, the replica should be covered with another sheet of glass or 
plastic to protect it from the dust. It can then be used as a lantern slide or for the 
preparation of photomicrographs. 

Figure 1 illustrates typical hailstones collected from two damaging storms which 
occurred in the Great Plains area of America during the Summer of 1959, They were 
collected in the field and shipped to us at the Loomis School more than a thousand 
miles away. They arrived 48 hours later in excellent condition. 

The first shipment was received by air express from Sterling, Colorado, nearly 
2000 miles from Conecticut. About twenty pounds of dry ice surrounded the 
shipment which was packed in crumpled newspapers contained within a double 
box of corrugated cardboard. The hailstones arrived 48 hours later, dry and in 
perfect condition with about half of the dry ice still remaining. 

The second shipment came from Grand Island, Nebraska. Although a similar 
packing arrangement was suggested, the packer placed the stones on top of fifty 
pounds of dry ice and surrounded the stones with crumpled newspapers. Some 
surface melting had occurred on the stones farthest removed from the dry ice but 
not enough to cause serious damage. This near disaster showed the importance of 
placing dry ice both above and below the specimens. In this latter case, more than 
half of the dry ice remained in the package after two days. 

By use of this technique and the cooperation of farmers, meteorologists, and 
other cooperative parties, the damaging products of hailstorms become readily 
available to the atmospheric physicist for research studies. 

The composite photo in Figure 2 illustrates an assembly of replicas prepared in 
the manner described. Figure 3 illustrates the detailed structure observed in the 
vicinity of a zone of concentric banding. Figure 4a illustrates a replica photographed 
normally and 4b the same replica as seen between crossed polaroids. The zones 
transmitting light represented areas containing voids which may be the same as 
those found to contain liquid water by List. These voids develop with the aging 
of hailstones at temperatures between —15°C and —25°C and may represent 
areas containing impurities. These and other interesting structures and phenomena 
observed during the current studies will be discussed in more detail in a subsequent 
paper. 

The replica method for studying hailstones has a number of interesting and 
attractive features. It permits the preparation of large numbers of permanent 
cross sections of hailstones, it permits the ‘‘reconstruction’”’ of a hailstone since 
successive sections can be made from one side of the stone to the other, and because 
of the ease of preparation and the relative inexpensive nature of the equipment 
needed, the method should lead to increasing the availability for study of a much 
larger number of examples of these interesting hydrometeors. 
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Zusammenfassung 


Es wird eine Technik fiir die rasche Zubereitung dauerhafter Abgiisse (Replicas) 
der Struktur von Hagelsteinquerschnitten geschildert. Das Verfahren diirfte sich 
als besonders férderlich in Studien iiber die innere Struktur solcher Partikeln er- 
weisen. Alle Vorkehrungen ké6nnen in einem Raum normaler Temperatur vorge- 
nommen werden, sofern eine Kiihltruhe mit einer Betriebstemperatur bis —20°C 
zur Vertiigung steht. 

Die Hagelschlossen werden auf einem Block aus plastischem Schaum verankert, 
mit einer Sage in zwei Teile geschnitten und die Schnittflachen mit grobem und 
feinem Sandpapier geschliffen und schliesslich iiber Samt geglattet, fiir welche 
Vorkehrungen nicht mehr als einige wenige Minuten benétigt werden. Hierauf 
werden die Eisk6rper fiir etwa 15 min bei einer Temperatur von —20°C gehalten, 
oder so lange, bis sich auf den polierten Schnittflachen infolge Sublimation (Ver- 
dampfung) Atzfiguren ausgebildet haben. 

Die geadtzten Eisoberflachen werden mit einer 5%igen plastischen Lésung von 
einer Temperatur von —5°C iiberzogen und nach Abtropfen der iiberschiissigen 
Lésung alsdann auf —20°C gehalten, bis das Lésungsmittel verdunstet ist. Die ent- 
standenen plastischen Filme werden hierauf von der Eisoberflache abgestreift und 
zwischen zwei Glasplatten gebracht. Sie kénnen in dieser Form entweder photo- 
graphiert werden oder schliesslich als Diapositive Verwendung finden. Lichtbilder 
von natiirlichen Hagelschnittflachen und solche von Abgiissen illustrieren die 
Leistungsfahigkeit des Verfahrens. 


(Received: November 4, 1959.) 
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Measurements of Two-Dimensional Zero-Lift Interference 
of Ventilated Boundaries at Speeds up to Sonic 


By Duanvapa M. Rao, and MarHAGONDAPALLY A. Ramaswamy, Bangalore, 
India?) 


Introduction 


The successful development of ventilated wall test-sections has not only elimi- 
nated the basic difficulties involved in transonic testing, but has also made possible 
the use of much larger relative model sizes than was hitherto considered advisable, 
in the conventional type of tunnels, from the viewpoint of wall interference. This 
latter fact is of particular interest to workers in most research institutions whose 
limited resources permit only small-scale facilities in the field of experimental high- 
speed aerodynamics. 

While it is known that, for the normal model blockage values, most ventilated 
test-sections have small interference at subsonic speed, for larger models the bound- 
ary interference may have a peculiar distribution over the model length, depending 
upon the wall geometry, which can lead to a cumulative interference of appreciable 
magnitude particularly with Mach numbers approaching unity. The interference 
pattern for Mach numbers greater than the model critical Mach number lies outside 
the scope of the linear theory, which also does not account for any viscous effects. 
The boundary layer is of considerable thickness particularly on perforated walls, 
and further, is subject to complex pressure gradients resulting from the model flow 
field being projected on the boundary with increasing intensity as the Mach number 
approaches unity, and its effect in modifying the boundary conditions represented 
by the slotted and perforated walls remains as yet relatively unknown. 

The purpose of this note is to present some results of experiments carried out at 
the High-Speed Aerodynamics Laboratory of the Aeronautics Department, Indian 
Institute of Science (during March/April 1958), with typical slotted and perforated 
walls in two-dimensional zero-lift flow with the specific aim of observing the 
interference characteristics of such boundaries. A half-nozzle arrangement was used 
in order to work at the maximum possible Reynolds Number in the small tunnel 
available, and a large relative model size was chosen to make the boundary inter- 
ference, and the variations in it, identifiable. By comparing the data from the 
present tests with ‘standard’, i.e. interference-free data for the same model [1]*), 
the variation of the blockage interference at mid-chord with free stream Mach 
numbers up to unity, as well as the chord-wise distribution of interference, were 
determined for three slotted and two perforated walls. 


Symbols and Definitions 
M7 Tunnel free-stream Mach number, 
Mr F Free-stream Mach number for interference free case, 


1) Department of Aeronautical Engineering, Indian Institute of Science. 
2) Numbers in brackets refer to References, page 334. 
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V. Tunnel free stream velocity, 

Mean cross-flow velocity through hole, 

V, Mean vertical velocity component at the boundary, 

q Tunnel free-stream dynamic pressure, 

Ap Pressure-drop through the wall, 

K,, = (46/9) (Vi/Voo)» 

K = (Ap/q)/(V,/V,.), 

Blockage interference factor: AM,,/M,, 7 = (M57 — Moop)/Mao 7 (Numerator is the 
difference between the tunnel and the free-air upstream Mach numbers correspond- 
ing to the same local Mach number on the model at a given point). 

: Total open area of slots/perforations : 

Ped MENS Nozzle (length x width) ; 


Experimental Details 


Figure 1 shows the details of the nozzle and the model, fitted to the 3” x1” 
blow-down wind tunnel of the High-Speed Aerodynamics Laboratory. The venti- 


UPPER NOZZLE 
BLOCK PLENUM PRESSURE TAPPINGS 


FLOW at VENTILATED PLATE (0-75° DIVERGENCE ) 
—— < 
MODEL 
a 
LOWER NOZZLE BLOCK 
STATIC ORIFICES 

DRILLED IN 
ARALDITE 

b 


0-2 0-4 0-6 0:75 0:9 


Figure 1 
a Details of test section. 6 Details of model showing static pressure orifice locations 
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lated test boundaries were represented by 1/8” thick aluminium slotted or perforated 
plates fixed to the upper nozzle liner (having a 1” deep plenum chamber) at a 
divergence angle of 0-75 degrees, relative to the lower liner surface, to compensate 
for boundary layer growth. 

The model, a 2”-chord 6%-thickness ratio circular-arc bump, was shaped from 
a block of ‘Araldite’ cast into the lower wooden nozzle liner, carrying the pressure 
lead-out tubes. The test Reynolds number based on model chord was about 1 mil- 
lion at M,, 7 = 1:0. 

All pressures were measured on mercury manometers, and converted to Mach 
number on the basis of free-stream stagnation pressure, using isentropic compres- 
sible flow tables. The Mach numbers are believed to be correct to + 0-005. 

Before conducting tests with the model, empty-nozzle flow was observed for 
all the ventilated walls, and the longitudinal pressure distributions were found to 
be satisfactory in all cases up to a Mach number of unity. With empty nozzle, 
there was a close agreement between the mean nozzle static pressure and the 
plenum chamber pressure for all Mach numbers. 

With the model present, the free stream Mach number in the nozzle, MW hes 
in all the runs was obtained by observing the static pressure distribution upstream 
of the model along the centre line of a side wall, and using a mean value of the 
static pressure in a region of nearly uniform flow. Due to considerable outflow 
through the walls induced by the model, there was a large difference between the 
plenum pressure and the mean free stream static pressure in all cases (e.g. see 
Pigure 3D). 

Visual flow observation was carried out during the tests using a 4-feet focal 
length concave mirror Schlieren system, and a series of Schlieren photographs were 
obtained using 1/100 s exposure on Kodak Tri-X film. Typical flow patterns are 
shown in Figure 4. 

The geometric details of the ventilated walls tested are given in the following 
table: 


No. | Type | Open Area Geometry 
1 Slotted 6-25% 2 slots x 1/32” 
2 Slotted 12:50% 4 slots x 1/32” 
3 Slotted 9-40% 3 slots x 1/32" 
4 Perforated | 18-20% 1/8” dia. holes 
| 5 Perforated | 10-20% 3/32” dia. holes. 


Results and Discussion 
(1) Axial blockage distribution 


The axial distribution of blockage interference, which is a measure of the flow 
distortion due to the tunnel boundary causing the flow past the model to depart 
from any possible free-flight flow pattern, has been theoretically shown to be sym- 
metrical for symmetrical models in solid, open and slotted test-sections, whereas 
_ for the perforated boundaries an unsymmetrical distribution of interference exists 
at subsonic speeds [2]. As pointed out by MarpeR and Woop [3], in the 
application of the thin aerofoil theory to the problem of perforated wall inter- 
ference, the imposition of the Kutta condition at the trailing edges of the solid 
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‘lifting’ portions of the wall (the ‘lift’ appearing through the pressure drop across 
the wall) causes a distinctive orientation in the flow. This leads to an asymmetrical 
distribution of the boundary interference velocity potential even when the obstacle 
shape is symmetrical. In physical terms, there is an inflow from the plenum in the 
model region, and a corresponding outflow from the nozzle upstream of the model, 
leading to reduced local velocities on the front part of the model and an increase 
over the rear. 

Further, the theory shows‘2] that with an increasing Mach number, this 
asymmetry is gradually reduced while the magnitude of the interference at the 
model increases until, with the Mach number approaching unity, the interference 
distribution tends towards that of the open jet (i.e. symmetrical) while its magni- 
tude at the model tends to infinity. 

Using the data measured in the present tests, an attempt was made to investigate 
the above theoretical predictions, and some results are shown in Figure 2, where 
the chordwise distribution of the blockage interference factor for a slotted wall 
(No. 3) and a perforated wall (No. 4) are compared for four Mach numbers. 


——O--- PERFORATED 44 4,—xX— SLOTTED + 3 


= SA - 


Figure 2 
Comparative axial blockage distributions for slotted and perforated boundaries 


. 
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It may be pointed out at this stage that in half-model tests such as the present 
ones (particularly in the two dimensional case) there is a certain amount of flow 
distortion at the model due to the displacement thickness of the boundary layer of 
the reflection-plane wall. The flow distortion indicated by the curves in Figure 2 
is therefore the combined effect of the wall-interference and the boundary-layer 
influence. The latter aspect may be considered briefly before discussing the bound- 
ary interference. 

The boundary layer effect in the present case consists primarily of the elimination 
of the front stagnation point, reduction of pressures near the leading edge, and a 
modification of the trailing edge condition due to wake interference. The major 
middle portion of the model chord however remains relatively unaffected due to the 
rapid flow acceleration and consequent thinning of the boundary layer [4]. In the 
present tests the flow over the middle 60% of the chord only has been considered, 
and therefore the data shown in Figure 2 may be taken to be substantially correct 
quantitative representation of the wall boundary interference. In any case, the 
boundary layer distortion effect remains constant as the ventilated walls are 
changed, and on this account the curves of Figure 2 will at least be expected to 
display the correct comparative trends. 

As seen in this figure, for the four values of the free stream Mach numbers in 
increasing order, there is a continuous increase of the mean blockage (of the free-jet 
type in all cases). The greatest difference between the interference distributions for 
the slotted and perforated walls, however, is apparent at M,,, = 0-805. For 
Mach numbers below this value, the magnitude of the mean interference over the 
chord is not large enough to show any systematic differences for the two kinds of 
boundaries, and for higher Mach numbers, as predicted by theory, the asymmetry 
of the blockage distribution for the perforated wall diminished so that the two 
curves show similar (i.e. fore-art symmetrical) form. It may be noted here that the 
theoretical value of the Mach number for maximum anti-symmetry (i.e. when 
K y1 — a = 2:56) as calculated for a doublet [2] for this particular per- 
forated wall is 0-93, a figure which is too high however, for the linear theory 
to be valid; (actually, as is shown in Figure 3 and will be discussed in the next 
section, at this Mach number for which the theoretical mid-chord interference is 
zero, the blockage factor as measured at this location is about 0-1). 

The blockage interference distribution curves for the other perforated wall 
(No. 5, not shown) were also seen to display similar characteristics. 

It may be mentioned that, considering the errors to which the blockage factors 
as computed here are subject (being the small difference between two Mach num- 
bers which are themselves correct only to + 0-005), the curves of Figure 2 show 
remarkable regularity. This further supports the belief that these data are reasonably 
accurate. 


(2) Blockage variation with Mach number 


The variations of mid-chord value of blockage interference factor with free 
stream Mach number for the different boundaries are compared with theory in 
Figure 3. It is seen that for the slotted walls (all the three of which have nearly the 
same theoretical interference), the observed values lie close to the theory [5] up to 
M37 = 98. as 

For the perforated walls, however, there is seen to be a substantial difference 
between theory*) and experiment even at subcritical Mach numbers, which may 


3) Theoretical curves calculated from References [2], [5] and [6]. 


332 Kurze Mitteilungen — Brief Reports — Communications bréves ZAMP 


SLOTTED WALLS 
+ OPEN AREA 12-5]. 
x ‘4 ” 9-447. 
A ” ” 6:257 


PERFORATED WALLS 
O OPEN AREA 18:27. 


Figure 3 
Mid-chord blockage variation with Mach number. 


be explained in terms of the rather large viscous effects expected for the particular 
walls tested. It has been observed that for perforated plates with hole diameter of 
the same order as the plate thickness, as in the present case, the pressure-drop/ 
cross-flow relationship, which is the fundamental factor affecting the interference 
characteristics of the wall, is irregular and non-linear [2]. In particular the slope 
Ky, = (Ap/q)/(V;,/Voo) remains small and considerably less than the theoretical 
value over the range of cross mass flows of interest. The final effect is to reduce 
the interference of a boundary which theoretically gives a solid-wall type of blockage, 
and increase that of a boundary having a theoretical interference of the free-jet 
type. Thus, as seen in Figure 3B, although theory predicts a solid-wall type 
blockage for both the perforated walls, the interference actually observed is nearly 


2 
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zero up to M,,; = 0:75, beyond which it is of the free-jet type, rising almost 
linearly to a blockage factor value of about 0-15. Further, the two perforated walls 
appear to have almost identical interference up to sonic speed, whereas the theory 
indicates a larger interference for the wall with a smaller open area. 

In Figure 3C, all the mid-chord blockage interference versus Mach number 
results discussed above have been collected together, and here is apparent a rather 
remarkable quantitative similarity between the interference characteristics of the 
slotted and perforated walls investigated in the present experiments. This is 
probably the consequence of very similar pressure-drop characteristics of the 
different walls, although they differ greatly in geometric configuration. Some 
evidence of this may be seen in Figure 3 D, which shows the relation between the 
nozzle free-stream static pressure for the slotted No. 3 and perforated No. 4 walls, 
obtained with the model present in the nozzle. The plenum pressure is lower than 
the free-stream value due to the drop caused by the outflow through the wall 
induced by the model. This pressure drop is seen to be indentical for these two 
walls, which, out of the five walls tested, are the most different from each other in 
geometric configuration. 


(3) Schlieren flow observation 


Two of the Schlieren photographs obtained during the present tests are shown 
in Figure 4. These were selected mainly to show the following features: 


Perforated Wall No. 4 Slotted Wall No. 3 


Se ae 


7 7 


Figure 4 
Schlieren pictures of flow past 6% bump with slotted and perforated boundaries M,, 7 & 1:1. 


(a) For the perforated wall, the intensity of disturbances introduced by per- 
forations into the nozzle flow, ee 

(b) For the slotted wall, the shear layer in the plenum chainber originating at 
the point of contact of the bow wave with the wall, indicating a strong outflow 
at this point, and 

(c) The highly distorted bow wave in both cases. 


Conclusions 


(1) Typical slotted and perforated boundaries, with area ratio ranging between 
6% and 18%, are found to yield negligible blockage interference (blockage factor 
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not exceeding 0-01) up to the model critical Mach number, with a model of blockage 
area as large as 8% in the two-dimensional zero-lift case. 

(2) For slotted boundaries, the available theory for two-dimensional blockage 
interference appears adequate up to the model critical Mach number. For the 
perforated boundaries, however, the actual cross-flow/pressure-drop characteristics 
of the given boundaries must be used for calculating the interference, even at 
sub-critical Mach numbers. 

(3) The slotted and perforated walls tested show almost identical mid-chord 
blockage interference variation with Mach numbers up to unity, in spite of the 
differences in geometric configuration. Some evidence is indicated to suggest a 
correlation in terms of the similarity of the pressure drop characteristics of the 
particular slotted and perforated walls. 
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Zusammenfassung 


In dieser Arbeit werden einige experimentelle Resultate aus dem High-Speed 
Aerodynamics Laboratory, Bangalore, veréffentlicht. Untersucht wurde die Stérung 
durch die halbdurchlassigen Seitenwande des Windkanals, wie sie bekanntlich zur 
Verminderung der Wandinterferenz im transsonischen Gebiet zur Anwendung ge- 
langen. Als Testobjekt wurde ein halbes, symmetrisches Kreisbogenprofil verwendet. 


Der Einfluss von drei geschlitzten und zwei durchlécherten Wanden mit zwischen — 


6 und 18% offener Flache auf die Sperr-Mach-Zahl wurde im Gebiet M = 0,7 bis 
1,0 gemessen. Bei einer Mach-Zahl bis 1,0 haben sich die Stérungseigenschaften 
in der Mitte des Profils bei den fiinf untersuchten Grenzflachen trotz den Unter- 
schieden in geometrischer Konstruktion als beinahe identisch erwiesen. Die exper- 
mentellen Befunde werden mit den Ergebnissen der heute verfiigbaren Theorien 
verglichen ; fiir die geschlitzten Wande wird zureichende Ubereinstimmung bis zur 
kritischen Mach-Zahl fiir den Testgegenstand gefunden, wahrend fiir die perforierten 
Grenzflachen die Wirkungen der Reibungskrafte zu betrachtlichen Unstimmig- 
keiten fiir alle subsonischen Machzahlen fiihren. 


(Received: December 3, 1959.) 
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SSS ee a ae eee eee 
Varia — Miscellaneous — Divers 


ASICA 


Association Internationale pour le Calcul Analogique 
(Bruxelles, 50 Avenue Franklin D. Roosevelt) 


International Association for Analog Computation 
Internationale Vereinigung fiir Analogierechentechnik 


deren Prasident, Prof. J. HorFMANN, in Briissel ist, beschaftigt sich mit allen Fragen 
der Analogierechentechnik. Dazu gehéren zwei Untergruppen. Zur ersten zahlte 
der eigentliche Analogrechner, der fiir die Lésung sehr vieler und verschiedener 
Probleme verwendet werden kann und der aus einer grésseren Anzahl gleichartiger 
Elemente, wie Verstarker mit sehr hohem Verstaérkungsgrad, Integratoren, Schalt- 
elemente, wie Widerstande, Potentiometer, Kapazitaéten und Dioden, sowie aus 
nichtlinearen Elementen besteht. Diese werden entsprechend dem «Analogschalt- 
bild» zu Ketten zusammengefiigt und gestatten vorwiegend die Lésung von Pro- 
blemen der Schwingungslehre und der Regelungsautomatik. Fiir sinus- und nicht- 
sinusf6rmige Schwingungen, fiir Frequenzgangbestimmungen und fiir transiente 
Vorgange eignet sich der Analogrechner besonders gut. Seine Funktionsweise liegt 
dem Ingenieurdenken recht nahe. Die Resultate erscheinen in Form von Oszillo- 
grammen, bei welchen relativ leicht der Einfluss von Anderungen der Parameter 
beobachtet werden kann. Die Prazision des Resultates hangt entscheidend von der 
Genauigkeit der einzelnen Bauelemente ab, sie liegt in der Gegend von Prozenten 
oder Promillen. Der Preis der Analogrechner steigt enorm mit der gewiinschten 
Genauigkeit. Partielle Differentialgleichungen lassen sich nicht so bequem behan- 
deln. Der Vorteil der Digitalrechner ist deren hohe Prazision. 

Die zweite Gruppe von Analogiegeraten sind Modelle oder Simulatoren; sie 
werden fiir jedes einzelne Problem besonders hergestellt. Dazu gehéren der elektro- 
lytische Trog, die elastische Membran sowie zwei- und dreidimensionale Netzwerke, 
mit denen partielle Differentialgleichungen behandelt werden kénnen. Zur Lésung 
von Potentialproblemen, z. B. der Feldverteilung in einer Triode, sowie in einem 
Transistor oder der Stroémungslinien um einen Fliigel eignen sich elektrische Netz- 
werke (rheoelektrische Analogie). 

Die ASICA fiihrt alle 3 Jahre einen grossen internationalen Kongress durch, so 
1955 in Briissel, 1958 in Strassburg und im September 1961 in Belgrad, und gibt 
hieriiber ein Berichtswerk in franzésischer und englischer Sprache heraus. (Presses 
Académiques Européennes, Bruxelles). Daneben erscheint alle 3 Monate ein Heft: 


Annales del’ Association Internationale pour le Calcui Analogique, 
Presses Académiques Européennes, Bruxelles, 


in welchen Einzelprobleme der Analogierechentechnik, Informationen tiber die 
internationale Tatigkeit und die Literatur besprochen werden. 
Die Analogierechentechnik wird auch in der Schweiz und an der ETH in stark 
cmendam Masse Verbreitung finden. 
Ep. GERECKE, Vizeprasident der ASICA, 
Institut fiir Allgemeine Elektrotechnik der ETH. 
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Beitrage zur Physik und Chemie des 20. Jahrhunderts. Lise Meitner, Otto 
Hahn und Max von Laue zum 80. Geburtstag. Herausgegeben von O. R. FRISCH, 
F. A. Panetu, F. Laves und P. Rossaup (Verlag Friedr. Vieweg & Sohn, Braun- 
schweig 1959). 285 S., 89 Fig.; DM 29.50. 

Die vorliegende Sammlung von Aufsatzen ehrt die drei grossen Forscher Lise 
MEITNER, Otto HaHN und Max von Lave zu ihrem 80. Geburtstag. Es ist eine 
erstaunliche Tatsache, dass diese drei Jubilaen innerhalb eines Jahres stattfinden. 

Natiirlich ist es nicht méglich, auf dem zur Verfiigung stehenden Raum die 
25 einzelnen Beitrage auch nur aufzufihren. 

Erwahnt seien aber die Beitrage der Ziircher Professoren HEITLER tiber The 
Penetration of Gamma-Rays through Matter and the Development of Radiation Theory 
und Laves iiber Kvistallstruktur und Kristallchemie von Elementen und metallischen 
Verbindungen. 

Etwa ein Drittel der Abhandlungen ist historisch orientiert und beschwort jene, 
den meisten von uns unbekannte, Zeit herauf, da mit heute unglaublich einfach 
anmutenden Instrumentarien die experimentelle Grundlage der modernen Atom- 
physik geschaffen wurde. Die ungemein wichtigen Beitrage der drei Gefeierten zu 
diesen Grundlagen sind geniigend bekannt und brauchen hier nicht aufgefiihrt zu 
werden. Es geniigt, darauf hinzuweisen, dass ohne diese Leistungen kaum eine der 
iibrigen Abhandlungen hatte geschrieben werden kénnen. 

Das Spektrum dieser Arbeiten ist ausserordentlich breit und reicht von der Geo- 
chemie bis in die Biologie, etwa zur rontgenographischen Abklarung des Baus von 
komplizierten organischen Molekiilen. 

Mit der Redaktion des Bandes waren O. R. Friscu, F. A. PANETH, F. LAVEsS 
und P. RosBaup betraut. Sie haben mitgeholfen, ein eindriickliches Denkmal zur 
Physik und Chemie der ersten Halfte dieses Jahrhunderts zu schaffen. R. Jost 


Handbuch der Kaltetechnik, Bd. VII: Sorptions-Kdltemaschinen. Von 
W. NIEBERGALL. Herausgegeben von R. PLANK (Springer-Verlag, Berlin, Géttingen 
und Heidelberg 1959). 540 S., 258 Abb.; DM 96.-. 

Jeder, der sich mit der Berechnung, Konstruktion und Anwendung von Sorp- 
tionskaéltemaschinen befasst oder auch nur diese didaktisch so ergiebige Maschinen- 
gattung etwas ausfiihrlicher im Unterricht behandeln méchte, wird es freudig 
begriissen, dass nun aus der Feder eines langjahrigen Mitarbeiters der Firma Borsig 
und eines der besten Kenner dieses Gebietes eine umfassende Darstellung der ver- 
schiedenen Typen und Anwendungsmoglichkeiten vorliegt. 

Nach einer ausfiihrlichen Behandlung der theoretischen Grundlagen mit An- 
gaben tiber Energie und Kiihlwasserverbrauch werden in einem zweiten Teil die 
fiir kleine Kalteleistungen — besonders fiir Kiihlschranke — geeigneten Apparate 
besprochen. Im dritten Teil geht dann der Verfasser auf sein ureigenes Arbeitsgebiet, 
die Sorptionskaltemaschinen fiir grosse Leistungen, ein. 

Wertvoll sind vor allem viele Angaben iiber die in der Praxis iiblichen Werte 
der umlaufenden Mengen, der erforderlichen Austauschflachen, der tiblichen Driicke 
und Konzentrationen sowie auch das Eingehen auf Wirtschaftlichkeitsfragen, die 
Inbetriebsetzung und die Behebung von Stérungen. Diese waren im bisherigen 
Schrifttum teilweise iiberhaupt nicht, teilweise erst nach langem Suchen zu finden. 
Wenn auch an einigen ganz wenigen Stellen eine noch etwas pragnantere Ausdrucks- 
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weise zu wiinschen ware und zum Beispiel auch die Behauptung, dass «die Ausfiih- 
rungen fiir kleine und kleinste Kalteleistungen, wie sie fiir Haushalt und Klein- 
gewerbe mit meist periodischer, aber auch manchmal kontinuierlicher Betriebsweise 
(jedoch ohne Fliissigkeitspumpe) verwendet werden» (S. 3), nicht ganz zutreffen 
diirfte, so vermittelt das Buch doch einen vorziiglichen Einblick in dieses auch vom 
Standpunkt der Verfahrenstechnik so interessante Gebiet. 

Das fliissig geschriebene Buch darf gerne all denen empfohlen werden, die 
irgendwie mit der sorptiven Kalteerzeugung zu tun haben. P. GRASSMANN 


Rheology, Theory and Applications. Bd. I. Von F. R. Erricu (Academic 
Press Inc., New York 1956). 761 S.; $20.-. 

Rheologie bedeutet dem Sinne nach Fliesskunde und umschreibt eine For- 
schungsrichtung, die die Begriffe «plastische Formung» und «Fliessen zaher, insbe- 
sondere nichtnewtonscher Medien» umfasst. Die Arbeitsrichtung ging aus der 
Kolloidchemie hervor, wandte sich aber in der letzten Zeit stark zur Werkstoff- 
physik hin und wurde durch die Griindung rheologischer Gesellschaften in verschie- 
denen Landern ausserordentlich gefordert. Wahrend der elastische Verformungs- 
zustand einerseits wie auch das Stro6men newtonscher Medien anderseits einen 
gewissen theoretischen Abschluss erfahren haben, liegt zwischen diesen beiden 
Grenzfallen ein Forschungsfeld, das von einer umfassenden theoretischen und 
experimentellen Durchdringung noch weit entfernt ist. E1ricH hat sich die Aufgabe 
gestellt, in einem dreibandigen Sammelwerk bedeutende Monographien aus diesem 
Gebiet zusammenzustellen. Der erste Band umfasst unter anderem eine Darstellung 
der Beziehungen zwischen Spannung und Verformung bei der plastischen Form- 
anderung von Metallen (D. C. DRUCKER), erganzt durch ahnlich gerichtete Arbeiten 
(W. PRAGER, J. FLEEMANN, G. J. DiENEsS und R. B. Dow). Anschliessend werden 
Fragen behandelt, die der Kolloidphysik und dem Verhalten der Kunststoffe naher- 
stehen; eine sehr interessante Ubersicht iiber akustische Eigenschaften von Fliissig- 
keiten (R. B. Linpsay) beschliesst den ersten Band. 

Wenn auch dieses Werk in erster Linie fiir den engeren Kreis der Fachleute 
bestimmt ist, so vermégen doch die einzelnen Arbeiten der chemischen Technik 
und der Technologie der Kalt- und Warmverformung von Metallen praktisch aus- 


wertbare Hinweise zu vermitteln und das Interesse weiter Kreise zu begriinden. 
J. R. DE FRIES 


Théorie des circuits non-linéaires en régime alternatif. Von VITOLD 
BELEVITCH (Gauthier-Villars Editeurs, Paris 1959). 296 S., 160 Fig.; bFr. 360.-. ; 

Das vorliegende Buch bildet gewissermassen die Erganzung zum friher erschie- 
nenen Buch des gleichen Autors Théorie des circuits de télécommunications, in 
welchem nur die Theorie der linearen Wechselstromschaltungen mit konstanten 
Elementen dargestellt wurde. Entsprechend dem Anwendungsgebiet, beschrankt 
sich der Autor auf die Betrachtung der Modulatoren mit Dioden zur Frequenzver- 
schiebung, der Demodulation mit Gleichrichtern, der Frequenzvervielfachung, der 
Oszillatoren und der Frequenzteiler. Allen Betrachtungen liegt die gleiche Voraus- 
setzung zugrunde, dass die Nichtlinearitét durch eine frequenzunabhangige reelle 
Funktion zwischen zwei Grossen (zum Beispiel Spannung und Strom) mathematisch 
formuliert werden kann. Dadurch werden Differentialgleichungen vermieden, sofern 
Stréme und Spannungen zeitlich periodisch verlaufen. Zur Vereinfachung der Ab- 
leitungen wird ferner angenommen, dass das frequenzabhangige lineare Element der 
Schaltung stark selektiv sei, so dass nur wenige Frequenzkomponenten geniigen, 
um den eingeschwungenen Zustand zu beschreiben. 
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Das Kapitel I resumiert die technischen und theoretischen Probleme, die sich 
stellen. In den folgenden Kapiteln II bis VII werden Gleichrichterstromkreise 
behandelt mit der Hauptanwendung der Frequenzverschiebung, insbesondere die 
lineare Theorie der Modulatoren zwischen selektiven Abschliissen, Gleichrichtung 
zwischen selektiven Abschliissen, nichtlineare Produkte sowie durch Unsym- 
metrien in Modulatoren erzeugte und eine allgemeine nichtlineare Theorie der 
Modulatoren. Das Kapitel VIII befasst sich mit der Frequenzvervielfachung zwi- 
schen selektiven Abschliissen. 

Die restlichen Kapitel (IX bis XIV) sind den Oszillatoren mit nichtlinearen 
Kennlinien gewidmet, und zwar werden behandelt: die freie Schwingung; spezielle 
Probleme, wo 2 Schwingungsfrequenzen méglich sind, die Synchronisation, die 
Frequenzteilung, Oszillatoren mit selbsttatiger Begrenzung durch Gitterstrom und 
eine allgemeine Theorie der Oszillatoren. 

Niitzlich sind die erganzenden Bemerkungen mit Hinweisen auf die Literatur, 
die am Schluss zusammengestellt zu finden ist. 

Die Behandlungsweise der Probleme kommt den Bediirfnissen der Praxis ent- 
gegen. Allerdings konnte dem Verhalten der Modulatoren bei hohen Frequenzen 
infolge Vernachlassigung der Sperrkapazitaten der Halbleiterdioden nicht Rechnung 
getragen werden. Aber auch ohne dies verdient das Buch volle Beachtung, besonders 
fiir den Unterricht an Hochschulen und fiir das Selbststudium, schliesst es doch eine 
spiirbare Liicke in der Fachbiicherei. H. WEBER 


Fonctions sphériques de Legendre et fonctions sphéroidales. Bande 1-3. 
Von Louis Rosin (Gauthier-Villars, Paris 1957, 1958, 1959). 209, 384, 289S., 
44 Fig.; fFr. 4000, 5000, 5500. 

Unter den besonders fiir die angewandte Mathematik wichtigen speziellen — 
Funktionen spielen die Legendreschen Kugelfunktionen und ihre Verallgemeine- — 
rungen als Lésungen der Schwingungsdifferentialgleichung eine hervorragende 
Rolle. Ihre in zahlreichen, weitverstreuten Publikationen untersuchten Eigen- 
schaften zusammenzufassen und iibersichtlich darzustellen, bedeutet eine immense 
Arbeit. Rosrn hatte die bis heute in franzésischer Sprache fehlende, moderne Dar- — 
stellung mit grosser Sorgfalt besorgt (fiir die er den Prix Poncelet 1957 der Académie 
des Sciences erhielt) und dabei nicht nur Bekanntes zusammengetragen, sondern 
auch etliche neue Beziehungen hinzugeftigt. Das Werk ist aber nicht einfach eine 
Formelsammlung, sondern eine zusammenhangende Darstellung, in der die Beweise, © 
sofern sie nicht zu umfangreich sind, ausgefiihrt werden; in den iibrigen Fallen wird 
dagegen auf die Literatur verwiesen. 

Um eine Idee vom Aufbau und Inhalt der drei Bande zu geben, seien hier 
wenigstens die Kapiteliiberschriften angegeben: I. Imtyoduction del’ équation différen- 
tielle associée de Legendre en Physique mathématique. Les polynémes de Legendre et les 
fonctions de seconde espéce, d’indice entier, positif ou nul. Il. Les fonctions associées 
de Legendre, Pr lt) et Qn” (mu), Pour n et m entiers, n= 0.111. Les harmoniques sphériques — 
de degré entier positif, négatif ou nul. IV. Les fonctions associées de Legendre de type 
général .V.Valeurs approchées des fonctions associées de Legendre : développements asymp- — 
totiques de ces fonctions et inégalités simples qu elles vérifient. VI. Séries de polynémes de — 
Legendre, de fonctions associées de Legendre et de polynémes de Laplace. VII. Formules — 
et théorémes d’addition pour les fonctions de Legendre. VIII. Les zéros ou vacines des 
fonctions Pi" (u) et QO” (u). 1X. Applications des fonctions de Legendre aux surfaces de 
vévolution, avec des systémes de coordonnées curvilignes orthogonales, autres que les 
coordonnées sphériques. X. Fonctions qui se vattachent étroitement a celles de Legendre: 
polynomes et fonctions de Gegenbauer ou ultrasphériques. Fonctions sphévoidales. — Auf 
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die drei Bande verteilen sich respektive die Kapitel I bis III, IV bis VI und VII bis 
X. Die einzelnen Kapitel enthalten auch zahlreiche bibliographische Hinweise. 
Die Gegenbauerschen Funktionen werden weniger ausfiihrlich behandelt, indem 
nur ihre wichtigsten Eigenschaften angefiihrt werden, und im Falle der Spharoid- 
funktionen wird nur eine Ubersicht (von etwa 50 Seiten) der Hauptergebnisse 
geboten, im tibrigen aber auf das bekannte Werk von J. MerxNner und F. W. 
SCHAFKE, Mathieusche Funktionen und Sphéroidfunktionen (Berlin 1954), verwiesen. 
Der erste Band enthalt im Anhang je die ersten Legendreschen Polynome und die 
zugeordneten Funktionen erster und zweiter Art; der dritte ein bis 1958 reichendes 
Verzeichnis numerischer Tafeln. — Fiir den praktischen Gebrauch wird man jeden- 
falls auch das alphabetische Inhaltsverzeichnis (das man sonst in den franzésischen 
Werken meist vermisst) zu schatzen wissen. Fiir jeden, der mit den Legendreschen 
Kugelfunktionen zu tun hat, wird das Werk von Rosin ohne Zweifel von grossem 
Wert sein, findet er doch recht iibersichtlich alle wesentlichen Eigenschaften und 
Beziehungen zusammengestellt. E. RotH-DESMEULES 


Probleme des Kreuzstrom-Warmeaustauschers. Von H. Kiinr (Sprin- 
ger-Verlag Berlin, 1959). 83 S., 34 Abb.; DM 10.50. 

Das Problem der Berechnung von Kreuzstrom-Warmeaustauschern, insbeson- 
dere der Ermittlung des mittleren Temperatursprunges in solchen Apparaten war 
schon seit langer Zeit immer wieder der Gegenstand theoretischer Untersuchungen. 
Bekanntlich bereitet es nicht geringe mathematische Schwierigkeiten. Die grosse 
Anzahl der technisch wichtigen Schaltungsméglichkeiten und Varianten und die 
teilweise bestehende Unsicherheit iiber die zu treffenden Annahmen haben dazu 
gefiihrt, dass das Problem immer wieder erneut aufgegriffen worden ist. Das Buch 
von KUtut liefert einen neuen wertvollen Beitrag zu diesem Problemkomplex. 
Ktut konzentriert sich vor allem darauf, méglichst einfache, direkt verwertbare 
Angaben iiber den Warmeriickgewinn zu erhalten. Darunter versteht man das 
Verhdltnis der tatsachlich iibertragenen Warmemenge zu derjenigen Warmemenge, 
die unter sonst gleichen Verhaltnissen in einem idealen Gegenstromapparat tiber- 
tragen werden kénnte. Besonderes Augenmerk richtet der Verfasser auch auf den 
Vermischungseinfluss, welcher fiir die in solchen Apparatesystemen zu erwartenden 
Temperaturverhiltnisse von wesentlicher Bedeutung sein kann. Es werden Nahe- 
rungsverfahren angegeben, deren Brauchbarkeit sehr sorgfaltig tiberpriift wird, und 
zur Erleichterung der Rechnung ist eine grosse Anzahl von Kurventafeln der 
Arbeit beigefiigt. 

Das nur 83 Seiten umfassende Bandchen ist damit das Ergebnis einer viel 
umfangreicheren Arbeit, als man vielleicht zundchst vermuten wiirde. Dem Inge- 
nieur, der solche Warmeaustauschapparate auszulegen hat, wird daher das vorlie- 
gende Buch, welches bei aller theoretischen Strenge auf die Bediirfnisse der Praxis 
zugeschnitten ist, dusserst wertvolle Dienste leisten. W. TRAUPEL 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLticce und E. Creutz, Band 38, 2. Teil: Neutronen und verwandte Gammastrahl- 
probleme (Springer-Verlag, Berlin 1959). 843 S., 338 Fig.; DM 178.-. 

Dieser neue Band der Fliiggeschen Handbuchreihe enthalt nur zwei Artikel: 
E. Amatp1, The Production and Slowing Down of Neutrons, und U. Fano, L. V. 
SPENCER und M. J. BERGER, Penetration and Diffusion of X-Rays. Die Arbeit von 
Amatpt (659 S., 281 Figuren) zeichnet sich durch ausserordentliche Griindlichkeit 
aus. Ausgehend von historischen Bemerkungen iiber die Entdeckung des Neutrons, 
wird dem Leser in einem einfiihrenden Abschnitt die Bedeutung des Neutrons als 
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Fundamentalteilchen und als Kernbaustein dargelegt. Die Diskussion der Wechsel- 
wirkung von Neutronen mit Materie, das heisst von Neutronen ausgeloste Kern- 


reaktionen, elastische und inelastische Streuungen, die Erzeugung kiinstlicher — 


Radioaktivitaten usw., fiihrt schliesslich zur Besprechung derjenigen Elementar- 
prozesse, die fiir die Verlangsamung von Neutronen verantwortlich sind. Die 
Theorie der Abbremsung von Neutronen wird — entsprechend den durchaus ver- 
schiedenartigen physikalischen und mathematischen Gesichtspunkten — in zwei 
separaten Abschnitten behandelt: Neutronenenergien grosser als einige Elektronen- 
volt und Neutronen in der chemischen Region unterhalb eines Elektronenvolts. Ein 
weiterer Abschnitt orientiert iiber die Diffusion thermischer Neutronen, wobei neben 
der theoretischen Behandlung des Diffusionsvorganges auch der messtechnische 
Aspekt gebiihrend beriicksichtigt wird. Das Kapitel ber die Produktion von Neu- 


tronen beschrankt sich auf das Gebiet der Neutronenenergien von einigen Mega- — 
elektronenvolt. Die Erzeugung von Neutronen in Reaktoren sowie von Hoch-_ 


energieneutronen wird in anderen Banden dieser Handbuchreihe besprochen. Der 
besondere Wert dieses Kapitels diirfte in der ausserordentlich griindlichen Re- 
prasentation des zurzeit vorliegenden numerischen Materials in Form von An- 
regungskurven, Winkel- und Energieverteilungen usw. liegen. 


Der zweite Artikel von Fano ef al. (183 S., 57 Figuren) befasst sich mit der © 


Durchdringung und Diffusion von Réntgenstrahlen unter folgenden einschranken- 
den Bedingungen: Die Materie wird als amorphe Anhaufung von Atomen betrach- 


tet. Das diskutierte Energiegebiet ist auf den Bereich beschrankt, wo die Penetra- — 


bilitat korpuskularer sekundarer sowie Bremsstrahlungseffekte noch keine dominie- 


rende Rolle spielen. Diese Konzeption erscheint insofern zweckmassig, als vieles aus ~ 


diesem Kapitel ohne weiteres auf Neutronen iibertragen werden darf und anderseits 


die Durchdringung hochenergetischer Strahlung durch Materie (Kaskadentheorie) — 


in einem anderen Band von G. MoLi£ReE besprochen wird. Beide Artikel sind in ~ 


englischer Sprache verfasst und von einem umfangreichen Literaturverzeichnis 
begleitet. F, HEINRICH 


Kinfihrung in den Warme- und Stoffaustausch. Von Ernst ECKERT 
(Springer-Verlag, Berlin 1959). 235 S., 171 Abb.; DM 28.80. 

Das Buch von Eckert, dessen 1. Auflage 1949 erschienen ist, liegt nun in iiber- 
arbeiteter und erweiterter Form vor. Schon die 1. Auflage zeichnete sich dadurch 
aus, dass auf sehr knappem Raum eine griindliche Einfiihrung in das Wissensgebiet 
des Warme- und Stoffaustausches gegeben wurde. Wenn auch die 2. Auflage dem- 


gegeniiber eine nicht unwesentliche Erweiterung erfahren hat, so handelt es sich 


doch immer noch um eine dusserst straff zusammengefasste Darstellung dieses 
heute so umfangreichen Wissensgebietes. Dass dabei trotz der Beschrankung des 


{ 
2 


Raumes eine derartige Klarheit und wissenschaftliche Strenge gewahrt werden — 


konnte, stellt dem Verfasser ein hervorragendes Zeugnis aus. Das Buch behandelt 


die drei Grundgebiete der Warmeiibertragung, namlich Warmeleitung, Warmeiiber- | 


gang und Warmestrahlung und schliesslich noch den mit der Warmeiibertragung 
eng verwandten Stoffaustausch. Spezifisch moderne Entwicklungen sind dabei in 
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die neue Auflage mit aufgenommen worden, so die spezielle Behandlung von — 
Schmelz- und Gefriervorgangen und des Warmeiibergangs bei hohen Gasgeschwin- — 
digkeiten und bei fliissigen Metallen. Im Kapitel iiber Warmeleitung ist ferner eine ~ 


kurze und klare Darstellung des Relaxationsverfahrens zu finden, eine mathema- 
tische Methode, die auch an vielen anderen Orten zweckmassig zur Anwendung 


f 


gelangt. Gewisse stérende Druckfehler, die in der 1. Auflage noch vorhanden waren, — 


sind jetzt ebenfalls verschwunden. Man kann diesem ausgezeichneten Buch nur die 
weiteste Verbreitung wiinschen. W. TRAUPEL 
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Mechanik 


von Prof. Dr. HANS ZIEGLER, Professor an der Eidgenéssischen 
Technischen Hochschule in Ziirich. 


Band I: Statik der starren und fliissigen Kérper sowie Festigheitslehre. 


Dritte, neubearbeitete Auflage. 
246 Seiten mit 389 Figuren. 
«Lehr- und Handbiicher der Ingenieurwissenschaften», Band 5 


Inhalt: 1. Statik der starren Kérper. 2. Statik der Blnseigecelton: 
3. Festigkeitslehre. 


Das Buch stellt eine moderne Neufassung von Band I der 
Mechanik von E. Meissner und H. Ziegler dar, deren zweite 
Auflage vergriffen ist. Dabei ist vermehrtes Gewicht auf die 
analytischen Methoden im Vergleich zu den graphischen und im 
Kapitel tiber Festigkeitslehre auf das nichtelastische Verhalten 
der Baustoffe gelegt worden. Ferner sind jetzt den einzelnen ae 
Abschnitten Ubungsaufgaben beigefiigt. 


Statistische Methoden 
der Populationsgenetik 


Ein Grundriss fiir Genetiker, Agronomen und Biomathematiker. ay 


Von Prof. Dr. Henri Louis LE Roy, ausserordentlicher ‘eke 
Professor fiir Biometrik und Populationsgenetik an der : 
Eidgenéssischen Technischen Hochschule, Ziirich. ; 4 


391 Seiten mit 65 Figuren. 

Reihe der experimentellen Biologie, Band 15. 

Sammlung «Lehrbiicher und Menographien aus dem Gebiete es : 
exakten Wissenschaften». <a 


Aus dem Inhalt: Die Merkmalspragung — Die Gen- und 
Umwelteffekte als Korrelationsfaktoren — Die Varianzanalyse 
als Grundlage fiir die Schatzung der Gen- und Umweltwirkung — 
Die kiinstliche Selektion — Literaturverzeichnis (mehrere 
hundert Literaturangaben tiber theoretische und praktische 
Studien). 


Uber die rein technische Anleitung hinaus, wie statistische 
Methoden auf genetische Untersuchungen anzuwenden sind, 
werden die genetischen Hypothesen aufgezeigt, welche von Fall 
zu Fall in die verschiedenen statistischen Methoden eingebaut sind. 
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